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Abstract. Although there is general agreement that a removal of classical gravitational singularities is not only a crucial 
conceptual test of any approach to quantum gravity but also a prerequisite for any fundamental theory, the precise criteria 
for non-singular behavior are often unclear or controversial. Often, only special types of singularities such as the curvature 
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singularities that arise according to the singularity theorems of general relativity. In these lectures we present an overview of 
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singularity theorems. This suggests possible routes for quantum gravity to evade the devastating conclusion of the theorems 
by different means, including modified dynamics or modified geometrical structures underlying quantum gravity. The latter 
is most clearly present in canonical quantizations which are discussed in more detail. Finally, the results are used to propose 
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semiclassical pictures of cosmological bounces. 
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1. OVERVIEW 

Physical theories are always idealizations without which the complexity of nature would be too great to fathom. 
Theoretical physics is, mostly very successfully, based on assumptions needed to formulate equations, find solutions 
and use them to describe, explain and further investigate physical phenomena. 

Sometimes, however, these assumptions may not be general enough for all purposes. When they are violated, the 
theory breaks down which mathematically appears as the development of singularities. An example is given by the use 
of continuous fields rather than discrete atomic structures in condensed matter physics. When fields vary too strongly 
on small length scales, such as in shock waves, singularities can occur in continuous field equations even though the 
basic, discrete physical description remains valid. Usually, deviations between solutions and observations increase 
before a mathematical singularity is reached. It is then clear that the approximate description can no longer be trusted 
beyond a certain point. But observations are not always available in such regimes where singularities are approached 
and an interpretation of mathematical singularities becomes more difficult. This is the case especially for gravity where 
observations of strong field regimes are lacking. 

Singularities in general relativity therefore play a special and dual role. First, the classical importance of singular- 
ities can be questioned since there are always assumptions behind special solutions or general theorems leading to 
singularities. But classical singularities in general relativity also provide an excellent chance to derive implications 
for the structure of space-time described by general relativity. When the theory breaks down, lessons for space-time 
structure result which can be especially important for the development of quantum gravity possibly replacing general 
relativity around classical singularities. 

We will first review the classical singularity theorems and sketch their main idea of proof. This will allow us to 
see which assumptions enter the theorems and what their main conclusions are. These theorems, rather than special 
singular solutions, define the singularity problem of general relativity. Their statements provide the measuring rod 
which any proposal for singularity resolution has to be compared with. 

The following section will deal with potential examples for singularity resolution which have been proposed in 
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FIGURE 1. Illustration of a smooth manifold dissolving into a classical singularity, or into a discrete space where quantum 
gravity takes over. 

quantum gravity, mostly in string theory and in canonical quantizations. The examples are not intended to be complete 
but to indicate the general types of ideas that have been put forward. (See also [JJ for a summary talk.) Here, quantum 
hyperbolicity will be formulated as a general principle and it will be shown to require characteristic properties of 
quantum gravity to be realized. 

Our specific formulation of the principle is worded in the language of canonical quantum gravity which is described 
in a subsequent section. We start with an explanation of the difficulties encountered in the first attempts of Wheeler- 
DeWitt quantizations, and show how their resolution naturally leads to loop quantum gravity. Quantum geometry and 
quantum dynamics in this framework are then discussed at length, as they provide the main pillars for any attempt to 
address the fate of classical singularities. 

In this framework, loop quantum cosmology has led to explicit constructions of dynamical laws from which non- 
singular behavior can be derived in several models. This is where the principle of quantum hyperbolicity is currently 
realized without counterexamples. Loop quantum cosmology thus provides the most general scheme of singularity 
removal available at present, and it can be used for explicit scenarios. 

In most cases, however, the basic description around a classical singularity requires deep quantum regimes which 
do not lend themselves easily to intuitive interpretations. It can thus be helpful to develop effective descriptions which 
capture some quantum effects but are otherwise based on classical concepts. This is available in semiclassical bounce 
pictures which provide examples of how singularities can be avoided through bounces in certain regimes. It also 
provides the basis for perturbation theory to compute phenomenological and potentially observable effects of metric 
modes and other fields traveling through a classical singularity. 

Although ideas in all five lectures are closely related, the text of any section can be read largely independently of 
the others. 



2. CLASSICAL SINGULARITIES 

General relativity^ describes the gravitational field as a consequence of the space-time structure determined in terms 
of the metric tensor as a solution to Einstein's equation Gab = ^^GTat sourced by the energy-momentum tensor Tab 
coupled through the gravitational constant G. In component form, these are coupled, non-linear partial differential 
equations of second order for the space-time metric gab- Unlike other fundamental field theories, they are generally 
covariant under arbitrary changes of space-time coordinates which shows that there is no background space and time 
on which fields are defined. This implies gauge symmetries for the space-time metric. 

Gauge theories in general are more systematically analyzed in a canonical formulation where one uses fields and 
their momenta rather than fields and their time derivatives. One is thus breaking up space-time tensors into spatial and 
time components which hides the underlying general covariance. To introduce this in general relativity, one foliates 
space-time into a family of spatial slices E, parameterized by an arbitrary time coordinate f . The space-time metric 
adapted to this foliation can be written as 

ds^ = -N^dt^ + qabidx" +N"dt){dx^ +N''dt) (1) 



^ We will mainly follow the notation used in f2l|. In particular, we adopt the abstract index notation where objects such as Tai, denote tensors 
rather than single components in specific coordinates. Expressions only true in certain coordinates are indicated with greek indices. We do not use 
different types of indices for space-time tensors and spatial tensors; it is rather clear from the context whether an object refers to a space-time or 
spatial manifold. Note also that N and N" denote different objects, a function in the first and a vector field in the second case, although traditionally 
the same letter is used. 



with the spatial metric qat on E,, the lapse function N and the spatial shift vector N". The components are all functions 
of spatial coordinates on E, as well as time t. When inserted in Einstein's equation, evolution equations in f , 



qab = ^^y^^ (^/'"ft ~ P'clah) + 2D^aNb) (2) 



and 



loitG Vdet<7 

result from the space-space components written in first order form for the spatial metric qab and its momentum 
p"'' = ^(K"^ - K^q"'' ) / 1 6 ttG which is related to extrinsic curvature Kab = (2N) ' ^ {qab - DaNt - DtNa ) of the spatial 
slices in space-time. In addition, there are constraints 



resulting from the time-time and time-space components of Einstein's equation in vacuum. In all equations. Da denotes 
the spatial covariant derivative compatible with qab, and ^^^Rah its Ricci tensor. It is clear from those equations that 
the components and N" do occur in the evolution equations of qab, but are themselves unrestricted (except by the 
condition that ([T]i must be a Lorentzian metric). 



2.1. Initial value problem 

We thus obtain an initial value problem only once lapse function and shift vector A^" have been specified 
throughout space-time as a gauge choice. (See [3] for a discussion of the types of initial value problems realized 
in general relativity.) Their equations can be interpreted as determining the manifold structure of space-time, which is 
most clearly seen when using the gauge source function F'' = S^^^vX' When conditions are imposed by fixing F'', a 

gauge is determined. Notice that the Christoffel symbol F^^ does not form a tensorial object, which we mark here by 
using the tilde on one index, and thus fixing its values even to zero restricts the choice of coordinate systems. 

Thus, fixing F'' implies conditions on the gauge functions and A^". Through the usual relation 
^vX ^ 18^'^ i^vgip + dxgvp - dpgyx ) between T^^ and gab, one has, for instance, 

d,N - N" daN = {K'' ~ n"Ta ) 

which is an identity between F'^ and the canonical metric components if F'' is kept free. (The vector field is the unit 
normal vector to the spatial slices.) If F'' is prescribed as a gauge choice, however, the equation becomes an evolution 
equation for A^ in this chosen gauge 
Moreover, we have 

- ^'''r;;,^; = -g^'-y^di = -v^v^x' (5) 

such that prescribing F'' poses "evolution" equations for space-time coordinates . A common choice is the harmonic 
gauge F'' = where coordinates are harmonic functions. Once coordinates on an initial spatial slice are chosen, space- 
time coordinates are determined as a solution of (|5]l by fixing the gauge source function. 

The space-time manifold as a topological set equipped with an atlas of coordinate charts is thus, to some degree, 
part of the solution as a consequence of general covariance. This is entirely different from other field theories for fields 
on a given background (metric) manifold as they are used for the remaining fundamental forces. This does not only 
give rise to complicated conceptual and technical issues when a quantization is attempted, as we will see later, but also 
to new physical features already present in the classical theory which are sometimes disturbing. We are able to derive 
properties of space and time themselves, and about their ends. Although not much is known about general solutions of 
Einstein's non-linear partial differential equations, there is one general feature common to most realistic solutions of 
general relativity; Space-time cannot be extended arbitrarily but develops boundaries where the classical theory breaks 
down. 



FIGURE 2. The ingoing normal family of a sphere is contracting because its cross-section area decreases along the affine 
parameter of the geodesic family. 

2.2. Singularities 

Since it is difficult to determine solutions or even asymptotic properties in general, a useful idea is to employ test 
particles as probes of possible space-time boundaries. One thus studies how test objects behave in a given solution 
to Einstein's equation and whether their motion, as described by the classical theory, has to stop at a certain point. If 
this occurs, the failure to move the test object further can only be attributed to a boundary of space-time itself since 
no interactions are included which could destroy the object. Quite surprisingly, this procedure allows far-reaching 
conclusions with only the slightest input from Einstein's equation flUt]. 

The precise criterion for a space-time singularity in this sense is geodesic incompleteness: space-time is singular 
if a geodesic, i.e. a word-line of a freely falling test object, exists which is not complete and not extendible. Thus, 
a curve e:.y ^ M exists whose tangent vector satisfies e" VqC^ = 0, i.e. it is a piece of an affinely parameterized 
geodesic defined on a proper subset C M which cannot be extended to be defined on a larger subset ^ D As 
the trajectories of freely falling particles, geodesies describe test objects subject only to the gravitational force. 

Before we sketch the main proof of a singularity theorem, we collect the typical assumptions and properties used: 

• Positive energy conditions, translated to positive curvature through Einstein's equation, imply focusing of families 
of geodesies and thus self-intersections and caustics. 

• Topological properties of space-time such as global hyperbolicity or spatial non-compactness then allow one to 
relate properties of geodesic families to space-time properties. 

• Appropriate initial configurations select physical situations, such as an everywhere expanding/contracting spatial 
slice (cosmology) or the existence of a trapped surface (black holes) where, with the preceding two assumptions, 
singularities are bound to occur 

The proof we sketch here demonstrates the importance of all these assumptions and illustrates how they could 
potentially be circumvented. To be specific, we focus on singularities as they occur in black holes, requiring the 
existence of a trapped surface as an initial condition. More details and different types of theorems can be found in 



2.2.1. Trapped surfaces 

A trapped surface is a compact, 2-dimensional smooth space-like submanifold T C M such that the families of 
outgoing as well as ingoing future -pointing null normal geodesies are contracting |7]. A geodesic family is defined by 
specifying a transversal vector field on a submanifold which uniquely determines a family of geodesic curves through 
each point of the submanifold in a direction given by the vector field at that point. For a compact 2-surface, we can use 
the inward pointing and outward pointing null normals as those vector fields, defining the ingoing and outgoing null 
geodesic families. As shown in Fig.|2] the ingoing family of normal geodesies is usually contracting in the sense that 
its cross-section area decreases, but conversely one intuitively expects the outgoing family to be expanding. A trapped 
surface requires even the outgoing family to be contracting and thus occurs only under special circumstances. 

To define expansion and contraction of families of null geodesies formally, we use their tangent vector field k^'. 
Then, expansion is defined as := V^fc". (This is similar to non-relativistic fluid dynamics where the divergence of a 
velocity field gives the expansion of fluid volume elements.) Similarly, one can use the tensor Bah = '^bka to introduce 



FIGURE 3. Light cone in Minkowski space. 



shear (Tab as its symmetric, trace-free part and rotation (Oah as its anti-symmetric part (paying due attention to the fact 
that k" is null in the precise definition which we are not going to need here). These tensors are subject to evolution 
equations along the family, following from the geodesic equation: 

k'V.Bat = k'{VhV,ka+R,iJkd) - Vb{k'V,ka) - {Vbk'){V,ka)+Rcb/k'Ki = -B' bBac+RdJk'k^ (6) 

where we used the Leibniz rule, the geodesic equation and the relation (V^Vi — V/,Va)a)c = Rabc''(>^d, valid for any 
smooth dual vector field (Oc, which introduces the Riemann curvature tensor Rab/ ■ 

A simple example for the computation of expansion Q is given by the light cone in Minkowski space whose 
generators define a family of null geodesies. The null tangent vector field is 

and in Cartesian coordinates of Minkowski space we have B^v ~ <5v^|i simply in terms of partial derivatives. Taking 
the trace, one obtains Q = 2sgn(f)/r which diverges at the tip of the cone where the null geodesies intersect and form 
a caustic. 

The relation between self-intersections of geodesic families and the divergence of their expansion is general; Take 
a geodesic family which initially fills all of space-time, i.e. which emanates transversally from a 3-dimensional 
submanifold. We can thus define three independent vector fields Kf — {d /dOi)" on a family of geodesies parameterized 
by coordinates cr,- of the 3-dimensional submanifold in addition to the null tangent k" = {d/dt)" along geodesies. A 
caustic (or conjugate point) forms if fcf becomes degenerate as a 3 x 3 matrix since the infinitesimal separation between 
different geodesies then vanishes in at least one direction. The matrix fc" is related to expansion by 

K^Bt = 4^bk" = k'^VbK^' 

using the commutation of coordinate derivatives, [k, (f,]" = k^Vb^ — fcf V^A:" = 0. Then indeed, expansion 

= = ( (C- > ) Ife'' V, K-f = V, log I det Kf I 

diverges when Kf becomes degenerate^ det Kf = 0. 

From the geodesic deviation equation (|6]l we obtain the Raychaudhuri equation as its trace: 

B = k'V.d = g"''k'V,Bab - - 5 0' - OabO"' + (OabOi"" - Rabk"k' . (7) 

Thus, if a family of geodesies is non-rotatingH which means (Oab — and curvature is non-negative as a consequence 
of energy conditions (such as the null energy condition Tabk"k^ > for all null vectors k") and Einstein's equation 



We have used (det ) 'Vjdet = [k ');,¥;, Kf which follows easily from det Kf = ^e„i,^^ie'j^' KfKjK^Kf . 

For a geodesic congruence one can see, using the Frobenius theorem, that this is reaUzed whenever the congruence is orthogonal to a hypersurface. 



FIGURE 4. Causal diagram in spherical coordinates showing the ingoing null geodesic family of a 2-sphere centered at the origin 
of spherical coordinates. The family has a focal point at the origin (left vertical line). Any point beyond the the focal point can be 
reached by a time-like curve from the initial 2-sphere. 



then 9 always decreases along geodesies in the family. This is the focusing effect of gravity and plays a major role in 
deriving singularity theorems. 

Quantitatively, we have 9 < —j9^ such that after integration 0^ ' > ' + j (f — fo) starting from initial values 9o 
at f = fo- This is the place where initial values for geodesic families enter: For negative 9o, 9 must diverge after finite 
time no larger than f — fo = 2/ 1 0o | ■ Geodesies must intersect before that time and form a caustic (focal point). 

2.2.2. From caustics to singularities 

Caustics are singularities where a geodesic family ceases to define a smooth submanifold of space-time. But they 
are not physical singularities since space-time itself is usually well-defined at points where light rays intersect. To 
relate the occurrence of caustics to space-time singularities we need one more ingredient beyond the initial conditions 
and positive energy assumptions together with Einstein's equation akeady used. 

The required basic statement from differential geometry is the well-known property of geodesies as extremizing 
arc-length. Space-Uke geodesies extremize arc-length — ^ gabe"e''dt among curves between two given points and 
minimize it unless there is a focal point between the two points. Similarly, time-like geodesies extremize proper time 
Ze = Jg. \/—8ahe"e'''it between two points and maximize it unless there is a focal point between the two points. For a 
null geodesic such an extremization condition is not possible since the norm of a null vector vanishes and thus any null 
curve has zero length. But there is an analog to the statement that minimization of arc-length by space-like geodesies 
or maximization of proper time by time-like ones ceases beyond focal points: Any point on a null geodesic beyond a 
focal point can be reached by a time-like curve. In a causal diagram in spherical coordinates, this is depicted in Fig.|4] 

According to the Raychaudhuri equation, together with the usual assumptions of positive energy, every null normal 
geodesic family starting from a trapped surface must develop a focal point since the initial expansions are negative by 
definition of a trapped surface. This is always realized for ingoing geodesies, but can it be possible for outgoing ones? 
There is no difficulty if space is compact, as shown in the left part of Fig. |5] since there will be another coordinate 
center or a periodic identification encountered by the outgoing family. In fact, for a compact spatial manifold one 
cannot clearly distinguish between ingoing and outgoing null normals of a 2-dimensional space-like submanifold. But 
for non-compact spatial topology the outgoing null geodesies continue to go on forever and can never be caught up 
with by a time-like geodesic. 

There is thus a contradiction between the Raychaudhuri equation, in the context of positive curvature and given 
the existence of a trapped surface, and space-time topology for non-compact space. At this stage, finally, topological 
conditions are needed and we are able to translate caustics of geodesic families into space-time singularities. The 
contradiction can only be avoided by concluding that incomplete geodesies exists. The outgoing geodesic family does 
not develop a focal point, despite of the Raychaudhuri equation, because it can simply not be extended arbitrarily. 
When we formulated the contradiction there was the hidden assumption that the null geodesies in the family can be 
extended arbitrarily, i.e. we assumed them to be complete. The only way to avoid the contradiction is to conclude that 
incomplete null geodesies must exist: space-time is singular in the sense of geodesic incompleteness (see the right part 
ofFig.|5). 



FIGURE 5. Causal diagrams of a spatially compact manifold with two focal points of null geodesic families (left) and of the 
outgoing null geodesic family encountering a space-time singularity in the spatially non-compact case (right). 



2.2.3. Scheme of singularity theorems 

The proof sketched above illustrates the general assumptions and conclusions used in singularity theorems: 

1 . Initial conditions ensure the existence of geodesic families with negative expansion. Typical cases are trapped 
surfaces, implying black hole singularities, or spatial slices whose expansion or contraction is bounded away 
from zero everywhere, giving rise to cosmological singularities. 

2. Using positive energy conditions together with Einstein's equation in the Raychaudhuri equation implies focusing. 
With the initial conditions specified, caustics develop in finite time. 

3. A caustic in general is only a "singularity" of the geodesic family, not of space-time. A singularity theorem finally 
results together with topological assumptions rendering a caustic into an obstruction to geodesic completeness. 

An important consequence is that dynamics of the gravitational field is not used very specifically, but only to translate 
positive energy conditions into focusing. 



2.2.4. Example: Schwarzschild geometry 

In general, it can be difficult to identify all trapped surfaces in a given space-time, but spherical trapped surfaces 
in a spherically symmetric spacetime are simple to detect. (See, e.g., lH, for non-spherical trapped surfaces.) We 
use this here as an example to illustrate the relation between the regions where trapped surfaces occur and space-time 
singularities. With a line element 

d/ ^ -N{r,t)^dt^+R{r,t)^dr^ + r^dQ.'^ 

formulated in the "area radius" r, we consider submanifolds defined by r — ro being constant. This is usually a time- 
hke submanifold as illustrated in the left part of Fig.|6l While the ingoing null geodesic family moves toward smaller 
r < ro and is contracting, the outgoing null normal geodesic family moves to larger r > ro and is expanding. Such 
a 2-sphere obtained as the cross section of a time-like constant-r surface is thus untrapped. When the submanifold 
r = ro is space-like, by contrast, both the ingoing and outgoing null normal geodesic families move to smaller r (or 
larger depending on whether the surfaces are future or past trapped); see the right part of Fig.|6l Any cross section of 
a space-like constant-r surface is thus trapped. 

Thus, whenever n„ = (dr)a is time-like, g^^nant, = R^^ < 0, we have trapped surfaces. (If we use a line element 
with non-zero shift A^,- = grrN'^, the inequality reads g'''' = {R^N^ — N}) /N^R^ < 0.) For the Schwarzschild solution 
R^^ = 1 — 2M/r, and any sphere with ro < 2M is trapped. To the future of this region there must thus be a singularity, 
as it is drawn in the usual conformal diagram Fig.|7] 



2.3. General situation and Alternatives 

Singularity theorems demonstrate the stability under non-symmetric perturbations of singularities explicitly seen in 
symmetric solutions, such as the Schwarzschild or homogeneous cosmological solutions. Historically, this played an 
important role in accepting the importance of singularities for solutions of general relativity. But as in any mathematical 
theorem, assumptions are certainly necessary. The most general theorems, those with the weakest assumptions, only 
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FIGURE 6. Constant area radius surfaces in splierical symmetry and null normal geodesic families, showing whether the surfaces 
are trapped or not. 




FIGURE 7. Conformal diagram of the Schwarzschild solution with time-like and space-like submanifolds of constant area radius. 



show that one incomplete geodesic exists. Moreover, the general structure (such as curvature divergence) is not 
illuminated at all. The significance of singularity theorems can thus be questioned, and indeed non-singular (though 
maybe not fully realistic) solutions exist even if positive energy conditions are assumed liollllll . 

Nevertheless, no general conditions for non-singular solutions are known, and thus singularities cannot be ignored 
in general relativity. Generic solutions have boundaries which cannot be penetrated by geodesic observers provided 
that positive energy conditions hold. General relativity is thus incomplete as it does not show what happens at and 
beyond boundaries of its solutions. Extensions of the theory are necessary. The key to solving the singularity problem 
is not to find non-singular solutions but to provide sufficiently general conditions under which non-singular behavior 
would be guaranteed. This is not available in classical gravity even if no energy conditions are imposed at all. 

Rather than using geodesic completeness as a criterion for non-singular behavior, one can use alternative conditions 
for non-singular space-times. One that is more physically motivated is generalized hyperbolicity lfl2|] which states 
that all space-times allowing a well-posed initial value problem for standard matter fields are to be considered non- 
singular. This uses physical, potentially fundamental fields rather than test particles. It is more general since (conical 
space-time) solutions are known which are geodesically incomplete but satisfy generalized hyperbolicity ifisl fl4ll . 
However, no general results are available at present while the usual cosmological and black hole singularities certainly 
present counter-examples. 

Instead of changing criteria for singularities, modifying gravity might lead to better situations: examples include 
alternative degrees of freedom such as test strings rather than particles in string theory or properties of quantized 
space-time in quantum gravity as they will be discussed in the next section. 



3. BEYOND GENERAL RELATIVITY 

Singularity theorems demonstrate geodesic incompleteness under quite general assumptions, based mainly on differ- 
ential geometry (the geodesic deviation and Raychaudhuri equations) in combination with positive curvature. This 
implies focusing of geodesic families which, together with topological conditions, result in space-time singularities. 
But special non-singular solutions do exist classically, so the question is not if one can avoid the conclusions by 
evading assumptions. What is missing is a general mechanism by which one can conclude non-singular behavior in a 
sufficiently general class of physical situations. 

Violating energy conditions is an obvious candidate to evade the usual singularity theorems, but even this does 
not work generally enough for the types of matter we seem to need for fundamental physics. Singularity theorems 
have been proven, with different assumptions, for instance in the context of inflaton fields which violate positive 



energy conditions Jistl- Moreover, dropping positive energy assumptions does not necessarily improve the situation 
but usually makes it worse as the development of sudden future singularities in so-called phantom matter field models 
shows 111611 . 

Other than that, only modifications of gravity itself rather than matter can help. But also this will be subtle because 
not much of general relativistic dynamics is being used in singularity theorems. 



3.1. Facets of the singularity problem 

The singularity problem is a complicated issue to be addressed in a more general theory of gravity, extending general 
relativity in a well-defined form. For any explicit discussion, the main difficulties are: 

1 . No general classification of singularities is available and many different types exist. 

2. Singularities are not always accompanied by unbounded curvature as in the best known examples. It is thus not 
sufficient to address only unbounded curvature because this is not even shown in singularity theorems. 

3. In fact, not all singularities should be resolved iITtIi . Some are useful to rule out negative mass, e.g. of the 
Schwarzschild solution. Such singularities are typically time-like rather than space-like which is an important 
additional property not covered in singularity theorems. 

In isotropic models, one can often construct bounce pictures where the volume of any solution is bounded away from 
zero. This can be achieved by appropriate modifications of the dynamics which avoid unbounded curvature. But such 
models do not serve as a general mechanism. 

Curvature singularities, although not necessarily implied by the singularity theorems, are the best known types 
realized in cosmology and black hole physics. They can often be dealt with in special ways which more or less 
directly ensure bounded curvature. In particular, large curvature implies high energy regimes of field theories defined 
on a singular gravitational background. As usually, field theories are expected to receive strong quantum corrections in 
high energy regimes which, when even the gravitational field is at "high energies" in the sense of strong curvature, may 
also modify gravitational dynamics itself. General relativity would then only be obtained as an effective description, 
or the small curvature limit of a suitable extension valid also at high curvature. 

In quantum field theory and condensed matter physics, effective descriptions on small energy scales or large length 
scales can be defined by integrating out "massive" or short wave length degrees of freedom which will become 
important at high energies. At low energies those modes can safely be ignored, but they will become relevant when 
typical energies reach their mass or when the curvature radius approaches their length scale. Ignoring the additional 
degrees of freedom in such a situation implies deviations, or possibly a failure of the low energy description in the 
form of singularities. It is then only the effective description which fails and appears to be singular while the more 
fundamental theory can (and should) be non-singular. 

This is an interesting and partially successful picture for which many examples are available. Most of these examples 
are motivated by special known solutions (mainly Friedmann-Robertson- Walker ones) and devised with cosmological 
bounces in mind. So far, they are not general enough to extend to more complicated, especially inhomogeneous, 
solutions. Moreover, they do not address the true problem of gravitational singularities: not all singularities have large 
curvature, and diverging curvature is not the basic mechanism behind singularity theorems. By addressing exclusively 
unbounded curvature one appears to be treating a symptom rather than the cause. 



3.2. Example: string theory 

String theory fl^ provides an example for a theory whose dynamics reduces to that of general relativity for small 
curvature and low energies but differs at large energies. In particular, string theory can be quantized perturbatively 
which is not the case for general relativity without high energy corrections. Conceptual features are, however, quite 
different in string theory. For instance, while solutions of general relativity are not just a metric perturbation but 
rather the space-time manifold itself, as discussed before, string theory in its present version formulates gravitational 
excitations on a given metric background. This certainly has implications for how generally the singularity issue 
can be addressed, keeping in mind that singularities in general relativity are understood as boundaries of space-time 
arising through the dynamical laws that govern its own structure. If a background manifold is put in from the outset, a 
discussion and resolution of singularities at a general level becomes impossible. String theory does, however, provide 



valuable insights into how special singularities can be resolved by new degrees of freedom. This occurs mainly through 
a different viewpoint on test objects replacing geodesies, and through candidates for massive degrees of freedom not 
contained in general relativity. Many examples are discussed, e.g., in the recent review ifioll to which we refer for 
details and more complete references. 

In string theory, one uses 1 -dimensional strings or higher dimensional branes rather than pointlike particles as basic 
objects. World-volumes of test strings or branes then replace geodesies followed by point particles as the submanifolds 
whose incompleteness would signal a singularity in the spirit of general relativity's singularity theorems. This has been 
shown to change completeness results and lead, in this sense, to more regular behavior especially for conical or orbifold 
singularities (which are not space-like). It also provides examples for new degrees of freedom which, when taken into 
account rather than being integrated out, are necessary for regular behavior ll20ll . 

This, however, has been difficult to extend to dynamical space-times such as those displaying curvature singularities. 
General arguments why such space-times provide a qualitatively different challenge have been presented in |21], and 
several so far unsuccessful alternative attempts can be found in Il22ll23|] . The main difficulty is that string perturbation 
theory generally breaks down in such space-times. (Perturbation theory might be useful for the singularity issue in the 
context of tachyon condensation [24], but has so far been used only for null singularities where curvature does not 
diverge.) 

A second source of additional degrees of freedom which could become light close to a classical singularity and help 
to resolve it are string winding modes around topologically non-trivial components of space-time (or brane separation 
parameters): their mass m oc R^^ is proportional to the inverse radius (measured in the background metric the strings 
are propagating in) of the compact direction they are winding around. This is highly massive, and thus negligible in 
low energy effective theory, if extra dimensions are small, but can become relevant in high curvature regimes. It has 
been shown that winding modes can easily lead to bounces [25], but typically only of some directions. With small extra 
dimensions one is necessarily dealing with anisotropic geometries such that a bounce in one direction does not imply 
a spatial volume bounded away from zero. Usually, only the compact dimension bounces. Moreover, by design special 
topologies or configurations are required for winding modes to exist which spoils prospects for a general mechanism 
as gravitational singularities also occur in simply connected space-times. 

These constructions referred to strings propagating in non-evolving backgrounds. It is much more difficult to find 
analogous mechanisms in dynamical space-times as they arise in cosmology or in the interior of black holes close to 
their singularities. In such a context, not only technical difficulties arise but also applications of low-energy effective 
actions, which are mainly being used to study the effect of massive degrees of freedom, are not always general enough 
1^^. While low energy effective actions are well-suited to study the propagation and scattering of fields which are not 
highly excited out of their vacua, dynamical space-times in quantum gravity have to include a gravitational state far 
away from its vacuum. Then, more general effective equations are necessary which, requiring good knowledge of the 
quantum gravity state, are more complicated to derive. 

The main difficulty is that in all examples one is still using test objects in a background space-time, such as new 
fields provided by positions of branes. This does not include the dynamics of space-time itself. Moreover, strong 
back-reaction effects occur ll27i l23il showing that gravitational dynamics is very relevant and that a pure background 
treatment is insufficient. The most detailed scenario in this context has been formulated for the Schwarzschild-AdS 
singularity, rather than a dynamical cosmological one. From an analysis of correlation functions of the conformal 
boundary field theory one can conclude that bulk properties such as horizons and singularities only emerge in the 
classical limit but are not present in quantum gravity ll28ll . 

Tight arguments have been put forward which indicate that bounces found in homogeneous models are unlikely 
to extend to inhomogeneous situations [29]. Although they have been discussed there mainly in the context of the 
AdS/CFT correspondence, some of the arguments are general enough to caution against direct generalizations of 
homogeneous results in any context, not just in string theory. Specifically, an upside-down potential for field modes is 
argued to arise in a boundary field theory description of the cosmological situation. Field modes in this potential will 
reach infinity in a finite amount of time, corresponding to the classical singularity. Since all modes are independent, 
they behave differently even if they started out in a highly correlated manner from an isotropic initial configuration. 
Quantum gravity, by way of a self-adjoint extension of the field Hamiltonian which leads to reflecting boundary 
conditions at infinite values of the fields, could make the behavior non-singular. But if this happens, the field modes are 
unlikely to return to an highly correlated state as the initial one. Thus, while evolution of the quantum theory continues, 
it does not easily lead to a classical bounce back to a classical geometry. The main property of inhomogeneities 
used in this argument is the large number of fundamental degrees of freedom, which all need to be re-excited 
collectively in a special way for a smooth bouncing geometry to result. Independently of the specific quantization 
of inhomogeneities, this is much easier to achieve in homogeneous models with a small number of degrees of freedom 



than in inhomogeneous ones. 



3.3. Geometry 

Singularity theorems are mainly statements about differential geometry as they refer to properties of geodesies on 
a curved manifold. Einstein's equation is used only at one place, relating positive energy to positive curvature, which 
then implies focusing effects in the Raychaudhuri equation. Focusing through positive energy is used in the most 
common theorems, but is not the only reason since singularities easily arise with violated energy conditions. 

For a general solution of the singularity problem one should thus focus on geometry, not just on dynamics. Although 
both are intertwined in general relativity, geometry determines which type of dynamics is possible such that it can be 
viewed as more basic. Rather than modifying gravitational dynamics in a given geometrical setting, using a different 
geometry could be much more successful to avoid the far-reaching conclusions of singularity theorems. 

An alternative geometry is automatically provided by background independent quantizations, such as canonical 
quantizations. Such theories are not based on objects in a background space-time but they quantize full metric 
components as the non-perturbative dynamical objects. For instance, most versions employ wave functions supported, 
e.g., on the space of spatial metrics q^b which is the configuration space of canonical general relativity. Geometrical 
objects then become operators acting on these wave functions with properties generally very different from classical 
smooth geometry. A new quantum geometry underlying gravity then arises with an entirely new setting for the 
singularity issue. Despite of the difference to classical geometry, for any consistent quantum theory of gravity smooth 
space-times have to be approached as a classical limit far away from singularities. But large deviations from classical 
behavior can occur around classical singularities, possibly resulting in regular equations. Quantum geometry then, if 
this picture is successful, provides links between classical parts of space-time which would otherwise be interrupted 
by singularities. 

Canonical techniques treat space and time differently because time derivatives of fields are replaced by momenta 
but their spatial derivatives are retained. This is also true in relativistic theories where manifestly covariant theories 
are rewritten canonically in a form referring separately to space and time. Although this hides the covariance of field 
equations, constraints ensure that solutions still respect the equivalence principle. 

In general relativity, a canonical formulation is defined by introducing a foliation of space-time into a family of 
spatial slices in terms of a time function f such that the sUces are E,: f = const. Moreover one chooses a time evolution 
vector field t" such that f " V„f ~ 1 which determines how points on different spatial slices are identified (along integral 
curves of t") to result in spatial fields "evolving" in coordinate time t. In addition to the time evolution vector field 
there is the geometrically defined unit normal vector field n" to the spatial slices in space-time. This allows one to 
decompose the time evolution vector field t" = Nn" +N", as illustrated in Fig. |8] into a normal part, whose length is 
given by the lapse function A^, and a tangential part given by the spatial shift vector A^". One can then split off the time 
components from the inverse metric g"'' by defining the inverse spatial metric q"^ = g"'' + n"nf' such that q"''ni, = 0. 
Thus, the spatial metric indeed has only components tangential to spatial slices. Solving the relation for 

N 

and inverting it results in the space-time metric dTJ in the canonical form akeady used before. As is clear from 
the construction, lapse function and shift vector arise through the space-time foliation and define coordinate time 
evolution through the way different spatial slices are identified. They are thus related to the space-time gauge but 
are not dynamical fields. This confirms the realization that they are not subject to evolution equations as their time 
derivatives do not occur in Einstein's equation. 

Evolution equations do result for the spatial metric qah{t), interpreted as a time dependent field through its values 
on different E,. In a canonical formulation we have first order equations (|2|i and (O for qui, and its momentum, related 
to extrinsic curvature K^,h — ^{J^ftlab — 2£>(aA^ft)) of E; (with denoting the Lie derivative along t"). These phase 
space coordinates given by spatial metric components q^i, and their momenta are subject to constraint equations ^ 
implementing their dynamics, as we will discuss in more detail in the next section. For now, it suffices to know that, 
when quantized, wave functions on the space of metrics (if a metric representation is chosen) arise which are subject 



FIGURE 8. Decomposition of the time evolution vector field t" into parts normal and tangential to a spatial manifold. 



to differential or difference (diffQ) equations depending on the quantization scheme. 

There are then no test objects in a background space-time, but wave functions or gravitational observables are 
the fundamental dynamical entities. Differential geometry becomes applicable only in a classical approximation to 
describe space-time, and only in classical regions would geodesies be defined at all. Geodesic incompleteness is thus 
inapplicable as a criterion for singularities in quantum gravity. This suggests a natural answer to why incompleteness 
occurs so generally in general relativity: geodesies themselves are only valid as long as a differentiable classical 
geometry can be assumed. (In fact, there are cases of geodesic incompleteness where simply the metric ceases to be 
differentiable. These are typically the examples where generalized hyperbolicity as a criterion leads to non-singular 
space-times although they would be geodesically incomplete.) When quantum geometry becomes relevant, geodesies 
stop and have to be replaced by something more appropriate. Only classical geometry would end, but not quantum 
gravity with its own version of quantum geometry. Such a scenario looks promising but it has to be developed and 
verified in detail, requiring explicit candidates for background independent quantum gravity such as a canonical 
quantization. The question then remains: How do we address or even define singularities in such a context? 



3.4. Quantum hyperbolicity 

Classical singularities refer to properties of spacetime as a metric manifold. Abstractly, they are thus identified 
through properties of the metric tensor on certain submanifolds, although an explicit classification in general terms 
is lacking. In a canonical quantization, classical singularities must then correspond to properties of states, or of 
suitable observables in a Heisenberg picture. The usual example is how a wave function in the metric representation is 
supported on certain submanifolds of the space of spatial metrics which classically would imply a singularity. Before 
we can even use quantum dynamics we must be able to provide an unambiguous one-to-one correspondence between 
classical singularities and metric tensors. This could to some extent be done in terms of curvature invariants, but must 
be more general since not any singularity refers to curvature. Such a classification is not available in full generality 
and is thus the major difficulty in any general discussion of the singularity problem in quantum gravity, iiTespective of 
the specific quantization approach followed. 

Although a general classification is not available, a new criterion for gravitational singularities becomes possible: 
If a state can uniquely be extended across or around all submanifolds of classically singular configurations they do 
not pose obstructions to quantum evolution in any sense. If, however, a state such as a wave function on the space 
of metrics cannot be extended uniquely across a classically singular submanifold, there would still be a boundary 
to quantum evolution and quantum space-time would remain incomplete. This viewpoint is the natural extension of 
generalized hyperbolicity from matter fields on space-time to the fundamental object of quantum space-time itself. We 
therefore call it the principle of quantum hyperbolicity. 

A verification in concrete ways requires quantum geometry for the location of potential singularities and quantum 
dynamics to see if states can uniquely be extended across such places. The set-up in this form is general without 
specific reference to unbounded curvature. It can potentially deal with the general types of singularities implied by the 
theorems. 



' Whenever we write "diff." we mean that differential or difference equations can occur depending on the context (the quantization scheme used) 
but the statement does not depend on the type of equation. 



FIGURE 9. A sketch of some part of superspace with classically singular regions. They could be finite regions or entire sub- 
manifolds splitting superspace into several disconnected components. The latter are most dangerous because they prevent physical 
information to be extended uniquely across that boundary. In addition, singular regions may be sharply defined submanifolds or 
washed-out regions due to quantum uncertainty. A detailed classification is not available in full general relativity but can often be 
completed in models which assume symmetries. 



States, as solutions to the dynamical equations of quantum gravity, thus have to be extended uniquely across classical 
singularities. As a general principle, it entails several sub-issues necessary for its verification: 

• First, phase space locations of classical singularities in terms of metric or curvature components have to be 
unambiguously identified. 

• The structure of submanifolds in phase space corresponding to classical singularities according to the first point 
then determines how the classical space of metrics is separated into several disconnected components. In some 
sense to be specified, any state has to be extended uniquely between any two such regions. This gives meaning to 
an extension of states across classical singularities. 

• Before this extension can be shown to exist uniquely, representations of states or relational observables must be 
chosen. The formulation of an extension is necessarily representation dependent because, for instance, no back- 
ground coordinates are available in which classical metrics qab{t) would approach a singularity. Any extension 
is thus provided in phase space variables (internal time) which are monotonic around a classical singularity. This 
selects representations in which one observes an extension of a wave function from one part of its support to 
another, rather than a superposition of different branches. 

• Finally, the dynamical equations of quantum gravity will be used to verify that all states can be extended in the 
way envisioned in earlier steps. For a general answer to the singularity issue it is important to prove this extension 
for all allowed states rather than for restricted classes such as those given by one specific initial condition. 

For a specific proposal of how quanturn hyperbolicity could be realized, one can use the usual Wheeler-DeWitt type 
wave functions on the space of metrics 1301 13111 . Then, metric variables would be used to provide extensions across 
a classical singularity. The first point above, classifying classical singularities through the metric behavior, cannot be 
performed currently in full generality, but is often easily available in symmetric models. One thus studies mini- or 
midisuperspaces of metrics respecting a certain symmetry and possible extensions of wave functions within this class 
of metrics. Also quantum dynamics can often be obtained explicitly in such models such that the program can be 
carried through. Also here, one has to be careful with interpreting results in any given model because the structure 
of classical singularities as well as quantum dynamics are truncated by restricting oneself to one class of metrics. 
But the investigation of several different models, covering different types of classical singularities, provides valuable 
information as to whether or not a general mechanism providing quantum hyperbolicity can exist. It also shows how 
details of the specific quantum theory of gravity used in the process matter which allows conclusions about general 
constructions. 

If realized, quantum hyperbolicity provides a quantum version of generalized hyperbolicity and thus deals with the 
well-posedness of evolution problems. If quantum equations of wave functions are well-posed, there are no boundaries 
for quantum evolution even where one would classically expect a singularity. Given that issues of hyperbolicity are 
very difficult in general relativity and generalized hyperbolicity at the classical level has not been studied much 
one could question the feasibility of quantum hyperbolicity as a verifiable criterion in quantum gravity. After all, 
quantum dynamics is expected to be much more complicated than classical dynamics from a mathematical point 
of view. In fact, quantum dynamics in background independent approaches has not even been fully formulated yet. 
Nevertheless, quantum theory allows an important simplification of testing well-posedness: If one uses a Schrodinger 
picture of states, dynamical equations are linear. This removes one of the major difficulties present in classical general 



relativity. The complicated issue is, of course, to extract observable information out of the states where the non- 
linearity of gravity enters. Moreover, solutions even to linear quantum equations can be difficult to find explicitly in 
inhomogeneous models with many independent degrees of freedom. But for hyperbolicity we only need to study the 
well-posedness of initial or boundary value problems without the need of having explicit solutions available. This is 
much more feasible for linear compared to non-linear equations. Then, given that states can be uniquely extended 
across classical singularities once quantum hyperbolicity has been verified, one is assured that observable information 
extracted from such extended states also reaches from one disconnected part of classical superspace to others. This 
is a hopeful sign that general statements about the singularity issue can be made in quantum gravity without having 
to face all the difficult issues, just as singularity theorems were derived in classical gravity without much knowledge 
about the general solution space. 

3.4.1. Examples: isotropic, homogeneous space-times and the BKL conjecture 

Nevertheless, simple examples demonstrate that there are many non-trivial issues in verifying quantum hyperbolic- 
ity. The general idea is best illustrated in isotropic models with a single metric component given by the scale factor 
a > 0, subject to the Friedmann equation (o/a)^ = ^^p{a,(p). Classical singularities for the usual matter contribu- 
tions, such as a scalar field, are reached at a = 0. This gives a simple and general identification of classical singularities 
for isotropic geometries. Since a is positive, the classically singular submanifold is a boundary rather than an interior 
submanifold. Wheeler-DeWitt wave functions in the metric representation take the form (p) with a scalar field (p 
assumed as the matter content. 

This provides the simplest situation, but is very special. Even keeping homogeneity, the situation changes consid- 
erably if anisotropy is allowed such as in the Bianchi I model. There are then three independent metric components 
fl/ > 0, / = 1,2,3, determining a metric ds^ — —dt^+Y,jaj{dx')^. Instead of the Friedmann equation, we have the 
constraint aid2 +0103 +020^ = which is solved by the Kasner solutions aj oc f"' where «/ are real numbers such 
that L/ OC/ = 1 = L/ 0!/ ■ These two equations have only solutions (except for the non-dynamical one corresponding to 
Minkowski space) satisfying — 1 < a/ < 1. One of them must then be negative while the other two are positive. The 
total volume is proportional to oiaifls °= t which vanishes at f = 0, corresponding to the classical singularity. Again, 
we have a simple characterization, but now the behavior of the metric components is very different: one metric com- 
ponent diverges at the singularity. The singular submanifold is thus not only a boundary but located at infinity. This 
makes investigating the well-posedness of initial value problems of wave function 1^(0 1,02, (23) in a neighborhood of 
the classical singularity more complicated. 

In any case, it demonstrates that conclusions drawn from isotropic models cannot easily be generalized because the 
structure of classical singularities themselves changes. On the other hand, there are indications related to the BKL 
conjecture |32] that anisotropic models are quite generic and provide crucial information even for inhomogeneous 
situations. In this context one is looking for generic asymptotic solutions close to a space-like curvature singularity at 
T ~Q, which are argued to be of the form 

ds^ = -dT^ + T^"! W (0)1 f + t2"2W + t2«3 W (0)3)2 . 

Thus, a homogeneous Bianchi model, whose invariant 1 -forms col, are used in the spatial part of the metric, is 
generalized by allowing position dependent Kasner exponents aj{x). For a Bianchi I model, for instance, one would 
simply have co' = dx' in terms of Cartesian coordinates x'. More generally, the 1 -forms co^ are left invariant 1 -forms 
on a 3-dimensional Lie group and thus satisfy the Maurer-Cartan relations do)' = — jCj^O)"' A 0)*^ with the structure 
constants Cj^ of the Lie algebra. 

Since inhomogeneities in this class of space-times are modeled by space dependent exponents o:/(x), there are 
gradient terms containing the spatial derivatives (9„a/ added to the homogeneous equations of motion. This can be 
derived from the Christoffel connection which, in addition to the terms occurring in a homogeneous model with 
constant a/, receives contributions 

where 5^„idx:"djc'^ = ^/(t^"''^' — T^"'){(i)')'^ is the difference spatial metric after introducing inhomogeneity. This 
term can easily be computed: 

- log T{d',daa^h) + ^adbOC^a) - 5''^5„,5rfa(„)) (8) 



which diverges logarithmically at T = 0. Thus, 3-curvature terms of ^^"^Rab/ resulting from this contribution, diverge 
as (logr)^, and the 3-Ricci scalar as T^"'(logT)^. 

This is to be compared to contributions of the homogeneous curvature scalar to see if inhomogeneities play a role 
for curvature and equations of motion. For a general Bianchi class A model with structure constants parameterized in 
the form Cj^ = ENmn^'^d'^ ||33|] we have the 3-Ricci scalar 

I-,-, 1 / n^flr n^a^ n'n^ n^Tv' ri^rr' 

^%hom = -T {^ + —^ + ^-2 2 2 

2 \q2q3 qm qm q^ qi qi 

for a diagonal metric of the form qu — q[i)5ij in the -basis. Assuming, without loss of generality, that tti < is the 
negative one of the Kasner exponents, the strongest divergence of '^'^hom is given by ^i/<72<?3 ~ -f2ai-2a2-2a3 Since 
2ai — 2a2 — 2a3 < 2a/ for all / if a\ is the negative exponent, the homogeneous contribution (|9]l is more divergent than 
the gradient contributions resulting from ([8]l. This gives rise to the conjecture that spatial derivatives are subdominant 
asymptotically close to the singularity and that a general solution behaves locally as a homogeneous model of the 
Bianchi IX form (n' = = = 1) where all diverging terms in '^'/?hom are present. 

However, the BKL conjecture, although by now strongly supported ll34l [35h . has not been proven, and it is even 
less secure that it will extend to a quantum version. Diverging curvature is relevant in the argument, and it thus refers 
explicitly to curvature singularities. Classical equations of motion have been used which themselves may change in 
such strong quantum regimes. Moreover, one must get close to the classical singularity in order to have subdominant 
spatial gradients. Bounce models, as they are often discussed, avoid this high curvature regime and cannot appeal to 
the BKL conjecture to justify a possible extension to inhomogeneous situations. Thus, inhomogeneous models as close 
to the full situation as possible must still be considered in detail. 




3.4.2. Boundaries 

As seen in the examples, when metric variables are used for a characterization, singularities usually appear at 
boundaries of superspace since dtiqab > by definition and singularities often have degenerate metric or triad 
components. This makes it difficult to extend states across such submanifolds; one could at most ask that the state 
at the classical singularity be well-defined and that, in a certain sense, the boundary does not pose an obstruction 
to physical evolution. This suggests, as it is often used, to impose boundary conditions on wave functions right at a 
classical singularity, with different motivations for specific proposals ifsil [36l[37i [ssll . If these are conditions which 
imply, for instance, wave packets being reflected off the boundary, one can argue that evolution is not interrupted. 
Such an interpretation is, however, not easy to justify in general. It is dependent on choices of the precise form of 
boundary conditions and difficult to make generic beyond isotropic models. Even if the wave function is not supported 
at a classical singularity, it does not mean that all extracted physical information is regular To test this one would have 
to construct observables and the physical inner product in order to compute expectation values. If all these quantities 
remain regular, one can conclude that the classical singularity has been resolved. But if only wave functions have 
been shown, or restricted by boundary conditions, to remain regular one is not in a position to decide about singularity 
resolution. For instance, while a wave function constrained to vanish at a = of an isotropic model taken together with 
the usual probabiUty interpretation (which is itself subject to difficult interpretational issue in quantum gravity) implies 
directly that any physical quantity is supported only away from zero, it does not lead to general lower bounds for a. 
While in any given state the expectation value of a in the usual inner product would have a positive lower bound, one 
can always choose a state where the expectation value comes arbitrarily close to zero. The boundedness of a is then 
only put in through selecting the particular boundary condition, but cannot be regarded as a consequence of quantum 
gravity. Moreover, the physical inner product in which one computes the expectation value might itself contribute 
a function of a diverging at a = 0. While one can sometimes compute the physical inner product in specific models, 
which does not only involve assuming symmetries but also selecting the matter content, there is no general information 
on its behavior Any interpretation of properties of wave functions which relies on the physical inner product is thus 
highly difficult to turn into a general argument to address the singularity issue. 

Probably the most advanced discussion in this context can be found in the application ll39l l40ll to gravitational 
collapse. It is shown there that a canonical quantization of the explicitly determined reduced and deparameterized 
phase space (i.e. an internal time has been introduced to describe physical evolution), which solves all the constraints, 
is non-singular in the sense that an initial wave packet not supported at the classical singularity always vanishes at 
the classical singularity. This provides a full analysis making use of Dirac observables and a physical inner product. 



It is by far the strongest result in this context. Nevertheless, the question of whether it can be extended to a general 
mechanism remains open. First, the direct use of Dirac observables, although necessary to discuss singularity removal 
possibly implied by boundary conditions of wave functions, makes it difficult to see what happens in other situations 
where Dirac observables are rarely known explicitly. The very fact that Dirac observables can be computed explicitly 
in this model may imply that it is quite special, including its dynamics. But more importantly, the mechanism in the 
collapse model depends on assumptions on a semiclassical initial state: One must assume that the initial state is not 
supported at all at the classical singularity but only at a peak value far away. Then, unitary evolution preserves the 
boundary condition which guarantees that the wave function will never be supported at the classical singularity. While 
this looks innocent, it makes singularity removal in this scheme unstable: If there is only a tiny contribution to the wave 
function not vanishing at the classical singularity, it would not change much of the initial semiclassical behavior. But 
this already spoils the preservation of the boundary condition and in general a wave packet approaching the classical 
singularity will not stay away from it completely. The mechanism, like others based on boundary conditions, is thus 
sensitive to precise details of physics at the Planck scale which have to be dealt with by extra assumptions. 

Classical singularities located at boundaries of the space of metrics thus do not offer an obvious general and 
verifiable way for regular behavior Quantum hyperbolicity does not require explicit observables but can be formulated 
directly for general states solving the dynamical constraints. Nevertheless, a location of classical singularities at 
boundaries seems to prevent the applicability of quantum hyperbolicity in any realistic sense because the choice of 
boundary data matters crucially. But one should note that the characterization of classical singularities as kinematical 
boundaries in this way depends on the variables used. The main reason for the location at boundaries was the restriction 
det<7„i > which obviously has to be imposed on the spatial metric. But geometry can just as well be described in 
triad variables, which is even necessary if fermionic fields are present. If a co-triad ej, is being used, defined as three 
co-vector fields e^, i — 1,2,3 such that Y^i^Wb ~ lab, the sign of detej, is relevant since it determines the orientation 
of space while det^„/, will still be non-negative. A surface dete^, = is then interior, not a boundary. 

A quick look at the models studied before reveals that this sometimes helps in rendering classical singularities 
interior submanifolds. In isotropy, there is again a single co-triad component e eM., e'^ — e5^, such that a = \e\. The 
classical singularity at e = is then an interior point. In these variables, classical singularities are boundaries not 
simply by definition of the basic variables but only for classical evolution which breaks down at e = 0. For quantum 
dynamics, the extendability of states thus becomes testable. 

In the Bianchi I model, however, we have e/ e M such that a/ = |e/|. While the two vanishing triad components of 
a Kasner solution would correspond to interior points, the diverging one implies that the classical singularity is still 
located at the infinite boundary of minisuperspace. A further reformulation alleviates this issue: we use a densitized 
triad Ef = |det(e^)|ef, with the triad e" inverse to ej,, e"ei, = 5/. Then, the single densitized triad component p eS. 
satisfies = |/7| in isotropic models where the classical singularity p = is still located in the interior In anisotropic 
models we have three components p' related to the metric components by oi = y/\p^p^ / p^ \ and cyclic. (There 
is a gauge transformation which changes the signs of any two components p' leaving the third fixed. Thus, only 
sgn{p^ p^ p^) — sgndet/if is physically distinguishable.) For a Kasner solution, p' °^ t^^"' where 1 — «/ > 0. All 
densitized triad components approach zero at an anisotropic classical singularity which is thus realized as an interior 



Typical homogeneous singularities occur as interior points of the densitized triad space. This is also true for the 
Schwarzschild singularity [42] and employing the BKL conjecture one can assume that general inhomogeneous 
singularities relevant for cosmology and black holes may have the same behavior. Given a candidate for quantum 
evolution, the extendability of wave functions can then be studied in finite neighborhoods. The choice of variables 
thus matters for quantum hyperbolicity which may seem an unwelcome dependence on coordinatization. But it is not 
only the characterization of classical singularities where the choice of variables matters but also, and even more so, 
for the success of a chosen quantization scheme. This determines which kind of quantum dynamics can be used with 
a set of states represented on the space of metrics or triads. It is thus an appealing possibility that the formulation of 
quantum dynamics may put restrictions on the choice of basic variables in a form which may or may not allow one to 
realize quantum hyperbolicity. The structure of classical singularities then becomes an important means to test general 
issues of dynamics in quantum gravity in a way nicely intertwining classical relativity with quantum geometry and 
dynamics. We will discuss these issues in detail in the following sections. 
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FIGURE 10. Isotropic minisuperspace in densitized triad variables p with a matter field (p. Quantum hyperbolicity requires a 
unique extension of wave functions across the line p = 0. 



3.4.3. Requirements on quantum dynamics 

The quantum hyperbolicity principle is testable as it requires crucial properties of quantum dynamics to be realized. 
Continuations of wave functions are relevant and thus the mathematical type of dynamical equations. On the other 
hand, the principle is insensitive to conceptual issues such as the interpretation of wave functions, observables or 
evolution which are largely unresolved in quantum gravity. Progress can thus be made on the singularity issue even 
before quantum gravity is fully understood. 

Using densitized triad variables in isotropic cosmology leads to a dynamical equation of the type A\j/p{(p) = Hxj/p {(p) 
with a diffll operator A on superspace and a matter Hamiltonian H which is usually a differential operator on the 
matter field also containing metric components. For quantum hyperbolicity we must be able, starting from suitable 
initial values at one side of /? = such as large positive p, to extend any solution uniquely across p ^0 as illustrated 
inFig.[IO] 

Any matter Hamiltonian H defined with a fundamental field rather than phenomenologically through an equation of 
state contains p^', such as in 

h=^-\p\-'^^pI + \p\'^^v{(p) 

for a homogeneous scalar with momentum pcp. The reason for this general behavior is that momenta in a relativistic 
canonical formulation caiTy a density weight. Since momenta appear in quadratic form in usual matter Hamiltonians, 
this requires an inverse of the determinant of the metric to get a well-defined, coordinate independent integration of 
the total Hamiltonian. This inverse determinant leads to in an isotropic reduction. For a well-defined evolution 

interpretation of the equation in p, the coefficients in the diff. equation must not diverge especially at = 0. This 
requires that a quantization taking into account gravity must lead to a matter Hamiltonian where the divergence of the 
kinetic term at /? = does not arise. One can interpret this as saying that quantum hyperbolicity in an isotropic model 
requires bounded curvature because classical curvature in an isotropic model behaves as an inverse power of the scale 
factor. This requirement is intuitively reasonable, but it is only secondary and derived from quantum hyperbolicity for 
isotropic models. 

This bounded isotropic curvature condition is much easier to check than quantum hyperbolicity and thus provides 
an easily accessible indication of its realizabihty. For the Wheeler-DeWitt equation (in an ordering as it occurs, for 
instance, in jillil ) 

\^i^iV\pW{p.<p)) = ^ {^\p\-''^^^^lf|^<p^-\p\'I^V{<p)^l,{p^<p)^ (10) 

with the ordinary quantization of p as a multiplication operator the condition is obviously not satisfied: the right hand 
side diverges generically at p = 0. Then, even with the classical singularity being interior in densitized triad variables, 
there is a break-down of the initial value problem ai p — Q. Wheeler-DeWitt quantum cosmology violates quantum 
hyperbolicity. 

Discrete approaches provide simple solutions to this problem if the continuous /7-space is replaced by a lattice not 
containing p ~0. Coefficients of the resulting difference equation are then never evaluated at p ~Q and no divergence 
arises. But as this avoidance of a singularity would merely be put in by hand by choosing a suitable lattice such a 
resolution would hardly seem satisfactory. (A more refined version of this possibility can be found in li45il46ll .') 



* Recall that this means differential or difference. 



FIGURE 11. Anisotropic minisuperspace with a classical trajectory and the singular hyperplanes. 



As we will see later, quantizations which pass the bounded isotropic curvature test of quantum hyperbolicity do 
exist. Bounded curvature in isotropic models thus provides a simpler but non-trivial test of quantum hyperbolicity. 
It is, however, not to be over-generalized, as it has occasionally been done in the recent literature, as anisotropic 
models show. In anisotropic models the classical singularity occurs at = ~ =0. This is a single point in the 
interior of 3-dimensional minisuperspace. But unlike the classical trajectory sketched in Fig.[TT] a wave function will 
be supported on 3-dimensional regions. Coefficients of the dynamical equations then can become singular even if a 
single p' vanishes. All three planes p' — are thus to be considered potentially singular for quantum dynamics. 

As before, coefficients of the diff. equation must not diverge for a well-posed initial value problem. Again, the danger 
comes from a matter Hamiltonian but possibly also from intrinsic curvature which can diverge in anisotropic models. 
Now, however, boundedness is not required on the whole minisuperspace; only a local version in any finite subset of 
the singular hyperplanes is necessary. Curvature can (and usually does) grow parallel to singular submanifolds, e.g. for 

^ while p^ and p^ remain small ll47ll48ll49ll . This shows that bounded curvature operators are not required in 
any non-isotropic model. An analysis of the singularity issue then has to refer to more detailed properties of quantum 
dynamics and the corresponding initial value problem. 

This is also the case in any inhomogeneous model. In such cases, many coupled diff. equations result, one for each 
spatial point as the equations are functional. Although equations remain linear, initial or boundary value problems 
are much more complicated and to be analyzed in detail. Without any symmetry it would not even be known how 
to locate classical singularities on superspace in a one-to-one manner. This complicates any scheme to resolve 
singularities, not just the principle of quantum hyperbolicity. Moreover, in full situations one may have to expect 
non-commutative metric operators to occur (see, e.g., ifsoll). If singularities are identified using metric components, 
such a non-commutativity would wash out classical results ifsTll . In such situations, only general principles can be 
applied which is reminiscent of the very general appearance of classical singularities in singularity theorems of general 
relativity. But even before such general situations can be addressed, in a combination of mathematical relativity with 
quantum dynamics, many non-trivial tests of quantum gravity arise from several symmetric or other models. The main 
criterion in those cases is whether or not the resulting initial value problem with given quantum recurrence relations is 
well-defined. 



4. CANONICAL QUANTIZATION 

To address the singularity issue from the point of view of quantum hyperbolicity, i.e. the unique extendability of 
states across classical singularities, a sufficiently detailed formulation of quantum dynamics is needed. Proper wave 
functions are required, independent of any background metric. Only then can quantum geometry be taken into account 
fully which on smaller scales or close to a classical singularity can be very different from classical geometry. Canonical 
quantizations allow systematic constructions of quantum field theories even without reference to a background metric. 
The metric itself can then be a full quantum operator and, e.g., fluctuate. Its quantum dynamics is relevant for the 
singularity issue. 

In canonical gravity, there is an infinite dimensional phase space, in ADM variables [52], of fields qat, the spatial 
metric, and momenta related to extrinsic curvature. The other components and A^" of the space-time metric ([T]) 
are not dynamical since and A'" do not occur in the action. Thus, according to the usual definition as the derivative of 
the action by time derivatives of fields, momenta p^ and pf^a vanish identically. Rather than determining the evolution 
of A^ and A^", their equations of motion = = — |^ and — pm — imply constraints on the phase space 



variables: the Hamiltonian constraint 



and the diffeomorphism constraint 

8H 

= — = IDhpl . 
5N" 

For a generally covariant theory the total Hamiltonian is a sum of constraints since no preferred time variable exists 
to which absolute evolution generated by a non-vanishing Hamiltonian would refer. For general relativity, we have the 
Hamiltonian 

H = H[N]+D[N"] = J d'xNix)^+ J A\N"{x)^^ (11) 

and the constraints determine the full dynamics. They constrain allowed values of the fields and their initial data, and 
they generate Hamiltonian equations of motion in coordinate time through Poisson brackets such as qab — {qab,H}. On 
the right hand side of those equations of motion lapse and shift occur through (fTTI) . For specific equations of motion, 
they thus have to be specified by choosing a space-time gauge. This determines which coordinate the dot in equations 
of motion refers to (e.g. proper or conformal time depending on whether = 1 or = a in isotropic cosmology). 

This is different in quantum gravity since there will be no reference to the space-time coordinates at all. Dynamics 
must be described in a gauge-independent manner rather than using space-time coordinates. While this is also possible, 
though complicated, to do in classical gravity, in background independent quantum gravity it is the only option. The 
tensors qab and p"'^ are then to be replaced by operators, acting on states such as \l/[qab] solving the infinitely many 
quantum constraints 

H[N]\l/ = D[N"]\l/ = for all and N" . 

The difficult part is to define precisely which function space "y/[qah]" refers to, how basic operators quantizing q^b 
and p"'^ are represented in a well-defined way, and how their products or even non-polynomial expressions in the 
constraints are being dealt with. 

All these steps simplify in homogeneous models such as in isotropic Wheeler-DeWitt cosmology. Classically, 
qab ~ cT'Sab just requires a single variable a to be quantized without any tensor transformation laws to be taken care 
of in this restricted class of coordinates respecting the symmetryQ Wave functions \j/{a) are simply square integrable 
functions of a single variable. The Wheeler-DeWitt equation is obtained by quantizing the Hamiltonian constraint: 

I d I d 

- ^ - -^a\j/°c //matter'/ 

a da a da 

where there is some freedom in choosing the operator ordering of a and d /da (the ordering here agrees with (fTOll). 
The diffeomorphism constraint D = vanishes identically. 



4.1. Index-free objects 

In general, a well-defined background independent quantization is much more complicated due to, for instance, 
non-trivial transformation properties of tensorial basic variables. A single component does not have coordinate 
independent meaning, but only the tensor q^b has. If we quantize single components, the question is which coordinate 
system an operator q^y should refer to. Properties of space are to be determined by states (e.g. through expectation 
values) only after operators have been defined. The space-time manifold in general relativity is part of the solution 
to Einstein's equation, whose coordinates follow from coordinates on an initial spatial slice once the gauge has been 
fixed fully. Thus, it becomes available only after the classical constraint and evolution equations have been solved. In 
quantum gravity, we have to turn the basic tensors given by the spatial metric and its momentum into operators before 
we can even formulate the constraint equations. There is thus no meaning whatsoever to an operator "q^v" because its 
classical analog q^v, when defined in one chosen coordinate system, will receive factors dx'^' /dx^ when transformed 



' In the spatially flat case the scale factor can be rescaled by a constant, which does however not introduce coordinate dependent pre-factors. 



to other coordinates. (These coordinates are only spatial because a canonical formulation deals with spatial tensors. 
Nevertheless, the classical tensors transform in this manner on any spatial slice, not just the "initial" one.) One can 
certainly avoid this by choosing a fixed set of spatial coordinates on any slice once and for all. But the resulting 
quantum theory would keep a trace of that choice and would be badly non-covariant. 

For this reason, no systematic quantization is known in ADM variables using the spatial metric and its momentum. 
No full quantum theory has been formulated in those variables but one has constructed exclusively models where 
one can reduce the metric variables to scalar quantities. Scalar quantities can then more easily be quantized since 
they do not receive coordinate dependent factors from the tensor transformation law. Examples include homogeneous 
models where the metric is determined by a finite number of spatial constants which do not transform under the 
allowed coordinate changes preserving the symmetry. (In fact, even such simple models can give rise to confusion 
from coordinate changes. An example is the scale factor of flat isotropic models which can be rescaled arbitrarily 
by a constant. This rescaling freedom cannot be appropriately dealt with in a quantization of the model unless the 
freedom is fixed from the outset.) Also in inhomogeneous models such as Einstein-Rosen or Gowdy models one can 
sometimes express the metric in terms of scalar fields on a given manifold f53','57l. But those are special properties not 
available in a similar form in full generality. It then remains unclear if those quantum models show typical properties 
of quantum gravity or special features used in their formulation. In particular, there is no well-defined relation between 
those symmetric quantum models and a potential full theory. 

To proceed, we thus need to reformulate unrestricted spatial geometries in terms of index-free objects. A successful 
classical reformulation starts by first removing one spatial index from <7„/, and (or K^b)- As already used, we 
introduce the co-triad ej, such that ej,e^ — qah- Here, / is just an index enumerating co-vectors and does not imply 
transformation properties under changes of coordinates. Moreover, the position of this index is irrelevant and will be 
summed over even when repeated in the same position. Similarly, we trade in an index / for a spatial index in extrinsic 
curvature by defining Kl, := e'-K^/h, using the triad, i.e. the inverse of ej,. We cannot proceed further in this way and 
remove all spatial indices because e"e', = 5'j would loose aU information about the metric. The triad, on the other hand, 

has all information about the metric. In fact, it has more freedom because an SO(3) rotation Rje[, leaves the metric qat, 
unchanged. This corresponds to the three additional components which the (non-symmetric) matrix e'^ has compared 
to the symmetric q^b- The new degrees of freedom are removed from the resulting field theory by imposing a further 
constraint on triad variables which has the form of the Gauss constraint in an SU(2) Yang-Mills theory. 

Using K'^, as one of the canonical variables leads to a momentum which is not exactly the triad but its densitized 
version, the densitized triad E" — |dete^|e". At this point we should recall that densitized triad variables were very 
convenient for the singularity issue because they implied positions of classical singularities in the interior of the space 
of geometries. At this point we see how the choice of basic variables also plays a role when defining a quantum theory. 
As we proceed, we will see that a reformulation of classical gravity in terms of index-free objects is possible precisely 
in terms of variables which use the densitized triad instead of the spatial metric or co-triad. 

To proceed with the definition of index-free objects, we will replace the tensor K^b by a connection. This will 
immediately suggest an index-free object, a holonomy, and has the additional advantage that spaces of connections are 
much better understood mathematically than spaces of extrinsic curvature tensors or metrics. Just as the metric defines 
a Christoffel connection, the co-triad defines the spin connection 

n, = -e'^-^.5(a[„4] + i.^.la[,4,). (12) 

As a functional of e[,, it has vanishing Poisson brackets with Ef and can thus not be used as a momentum replacing 
K'l. But the Ashtekar connection [55, A^, = F', — yK'^^ transforms as a connection, as any sum of a connection 
and a tensor does, and is canonically conjugate to the densitized triad because K'^ is. In the definition, 7 > is the 
Barbero-Immirzi parameter ll56ll57ll which does not have any effect in the classical theory but is important in quantum 
gravity. 

Loop quantum gravity is based on a canonical quantization of the phase space spanned by A'^ and Ef with 

{AUx),E^iy)} = 87zrG5^5}5ix,y) . (13) 

Its success relies on the fact that connections and densitized vector fields can easily be expressed in terms of index-free 
objects which can then be represented on a Hilbert space. Instead of connection components we use holonomies 



heiA) = ^exp / A^,T,e"df 



(14) 



for curves e with tangent vector e" and Pauli matrices T, = — jOy. As usually, ^ denotes path ordering along e of the 
non-commuting functions in the integrand. If holonomies are known for all curves e in space, the connection can be 
reproduced uniquely fs^. Similarly, we use fluxes 

Fs{E)^ jd^yn^Efx, (15) 

for 2-surfaces S with co-normal ~ j Sabc^"^ Itt ^ parameterization S:y i-^ x{y). Again, if fluxes are known for 
all surfaces S in space, the densitized triad is reproduced uniquely. 

Notice that no background metric is used in these definitions as the tangent vector as well as co-normal are defined 
intrinsically without reference to a metric. Moreover, there are no free spatial indices and the objects transform trivially 
under changes of coordinates. Instead, there is a representation of active spatial diffeomorphisms (p which move along 
the labels e i-^ (j){e) and S t-^ (j){S) of holonomies and fluxes. These objects can thus be represented on a Hilbert space 
without having to include coordinate factors in the tensor transformation law. 



4.2. Loop quantum gravity 

For a representation one also has to know the Poisson algebra of basic variables which is to be turned into a 
commutator algebra. Poisson relations of holonomies and fluxes define the holonomy-flux algebra in which no delta 
functions occur even though we are dealing with a field theory: integrating connections and densitized triads to obtain 
index-free objects has automatically introduced the correct kind of smearing. Any delta function present in the Poisson 
relation (fTsT i between A'^ and Ej is integrated out in the holonomy-flux algebra. A well-defined algebra results, ready 
to be represented for a quantization. This provides the kinematical setting of loop quantum gravity. 

There are then two routes for applications of this basic structure: it allows a formulation of quantum dynamics in a 
background independent manner which we are going to describe now. Moreover, the unique representation of the full 
theory induces representations of symmetric sectors derived from a full framework. This leads to symmetric models 
with a precise relation to the full theory, which was lacking in Wheeler-DeWitt quantizations. As we will see in the 
next section, the induced representation is inequivalent to the Wheeler-DeWitt representation even in the simplest 
models, with important implications for the singularity issue. 



4.2.1. Representation 

One can now construct a quantum representation of the smeared basic fields which arose naturally in providing 
variables suitable for a background independent quantization. After the basic objects have been represented, thus pro- 
viding a quantum theory, one can quantize and impose constraints as operators. This will ensure that only space-time 
covariant observables arise. From the form of the basic variables, an SU(2) connection and a densitized momentum 
field in Ashtekar variables, general relativity appears as a gauge theory subject to additional constraints. 

A representation is most easily constructed in the connection representation where states are functionals on the 
space of connections. This would not be available had we not been led to introduce connections rather than tensors in 
constructing index-free objects; the representation will thus carry characteristic features as traces of the background 
independent quantization. In a connection representation, holonomies act as multiplication operators. Starting from a 
basic state which, as a function of connections, is constant, we thus "create" non-trivial states which depend on the 
connection along the edges used in multiplicative holonomies. Although such states depend on the connection only 
along edges, the resulting states can be complicated with edges non-trivially being knotted and linked with each other. 
Moreover, edges can intersect each other giving rise to vertices in which more than two edges meet. The resulting 
space of states can be spanned by a basis of spin network states fs^ defined as 

reg eeg 

where g is an oriented graph collecting all the edges used for holonomies, with labels j^, indicating irreducible SU(2) 
representations p/^ in which edge holonomies are evaluated, and matrices Cv which ensure that matrix elements of 
holonomies are multiplied with each other in a way resulting in a gauge invariant complex valued function of the 




FIGURE 12. A colored graph representing a spin network state togetiier witii a 2-surface intersecting a link. 



connection. Since holonomies are multiplied with SU(2) group elements at the endpoints of their edges under a gauge 
transformation, the contraction matrices sit in vertices where different edges meet each other. 

A well-known example of a gauge invariant function of holonomies is the Wilson loop obtained by taking the 
trace of a holonomy around a closed loop. In the fundamental representation, we can reproduce this in the above 
language by introducing two vertices vi and V2 along the loop which split the loop into two non-selfintersecting 
edges ei and 62- We orient this graph such that both edges start in vi and end in V2- Their holonomies in the 
fundamental SU(2) representation are 2x2 matrices hj{A)^ which change under a local gauge transformation to 
g{vi)f)hj{A)'j^g^^{v2)^ where g{vj) S SU(2) are the values the local gauge transformation takes in the vertices. There 
are thus two factors in each vertex, such as g{vi)^^^g{vi)'^^ from the two holonomies, which in the final expression 
must cancel each other when contracted with the matrix Cy^ . Using the identity SABgc — (^^)a^bc satisfied for any 
SU(2)-matrix g, one can see that a matrix C(vi)/ib = Eab and similarly C(\'2)^^ = results in a gauge invariant 
function: C(vi)ABg(vi)^g(vi)£ = C{vi)cd- The resulting spin network state 

W,(A) ^e^ce^^. (/ii(A))^(/,2(A))g = ih,{A))i{h2{A)-'f^^tr{h,{A)h2{A)-') 

leads to the usual expression for the Wilson loop, where in the second step we used the same SU(2) identity as before. 
At intersection points of higher valence, one can use decompositions of tensor products of SU(2) representations into 
irreducible ones, i.e. the usual recoupling rules known from angular momentum in quantum mechanics, to find all 
vertex matrices C,, leading to gauge invariant results. 

Spin network states span the whole quantum representation space because the action of holonomies as basic 
configuration variables is complete. On any such state, holonomies are represented as multiplication operators as used 
in the construction of the states. Since these operators have to respect unitarity properties representing the classical 
reality conditions of A^,, an inner product of the representation space results. In this process spin network states turn 
out to be an orthonormal basis provided one chooses an orthonormal representation of the vertex matrices Q. (which 
for a given number of edges meeting the vertex and given edge labels form only a finite-dimensional space as it follows 
from the rules of tensor decomposition). When completed, this defines the kinematical Hilbert space of loop quantum 
gravity. 

Fluxes are conjugate to holonomies and thus become derivative operators on the Hilbert space. Using 

which is a consequence of the chain rule for A', derivatives acting on functions of holomies, shows that non-zero 
contributions result only if S intersects the edges of g. Indeed, the functional derivative of holonomies (assuming, to 
be specific, that the surface intersects the edge at its starting point) gives 

f 8h \ f (' 1 

Js'^^^""WJy) " jdtna{yy5[e{t),y)he = -lnt{S,e)x,he 

and thus a factor Int(5',e) of the oriented intersection number of the surface S and an edge e. Individual contribu- 
tions in the sum over intersection points are then determined by "angular momentum operators" (su(2) derivatives) 
(T,7i)^(9/(9/!g = Xx{Tihd / dh) acting on holonomies. Since these operators have discrete spectra and are summed over 
in an at most countable sum, all flux operators have discrete spectra. With fluxes being the basic operators representing 



spatial geometry through the densitized triad, discrete spatial geometry emerges from the construction without being 
put in in the first place. 

This representation is not only convenient to construct and to work with, it is also, under mild assumptions, the 
unique representation of the algebra of holonomies /z^. and fluxes on which the diffeomorphism group acts unitar- 
ily ll60[ I61II62I1 ■ Classically, diffeomorphisms act on holonomies and fluxes by moving the defining submanifolds, 
{he,Fs) I— *■ If this is required to carry over to the quantum theory, as it should be since any violation of 

unitarity of the diffeomorphism group would imply a breakdown of spatial background independence, no other repre- 
sentation is possible. As it happens often, the requirement of a symmetry reduces the class of available representations. 
With the large diffeomorphism group as a consequence of background independence, the representation appears to be 
selected uniquely. 



4.2.2. Quantum geometry 

While fluxes do not have direct intuitive implications for spatial geometry, they occur in more typical objects 
such as the area operator. The area of a surface S with co-normal n„ as used in the definition of fluxes is A{S) = 

Jgd^y^J EfriaEfriti. A quantization thus requires a product of flux operators which can be defined after regularization 
ll63i l6?l. Due to the square of triad components present in the classical expression, the quantum operator contains a 
square of angular momentum operators whose spectrum is well-known. This allows one to determine the area spectrum 
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valid for the case where no intersections between S and the graph occur in vertices of g. In the general case the 
spectrum is more complicated but also known explicitly. 

Similarly, a volume operator V is obtained by quantizing the classical expression V{R) = /^d^xy^detEj. Again 
after regularization, contributions now come only from vertices v E RD g whose values are constructed from the 
invariant matrices in the vertex ll63[l65ll . Although the spectrum is much more complicated to determine than the area 
spectrum and not known completely, it is discrete. 

The volume spectrum contains zero as a highly degenerate eigenvalue which is realized for instance, but not 
exclusively, in the case where no vertex lies in the region R. But also the total volume of the whole spatial slice 
has a highly degenerate zero eigenvalue even for vertices with arbitrarily high valence. Therefore, there is no densely 
defined inverse of V. However, when we come to matter Hamiltonians we will need the inverse determinant of E", 
such as in the case 



f , / I pl+E"Efda(pdi,(p I 
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of a scalar field. Since also metric components entering matter Hamiltonians have to be quantized in quantum gravity, 
we need an inverse volume operator. As the spatial volume vanishes usually at classical singularities, this issue is 
related to the singularity problem. In fact, in the context of quantum hyperbolicity we have already seen that such 
coefficients matter for the well-posedness of initial value problems of wave functions. 



4.2.3. Quantization and ambiguities 

Using 

\a'„ ( Vldet^ld^A - 27trGe'j''eabc^== (17) 

and approximating A^, by holonomies one can replace inverse powers of detE by positive powers in Poisson brackets 
ll66ll . Inserting holonomies and the volume operator, and replacing the Poisson bracket by {ifi)^^ times a commutator 
then results in well-defined operators with a classical limit as required. But the resulting operators are not identical 
to an inverse volume which does not exist in the quantum theory. Deviations between the classical inverse and the 
quantum behavior thus result which are most noticeable at small length scales. 
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FIGURE 13. Loop a used to quantize curvature components. 



With many different ways to re-write, e.g. using different representations for holonomies, giving the same classical 
expression but differing in quantum properties such as their spectra, the operators are subject to quantization ambigu- 
ities. As they are non-basic operators and in fact non-polynomial functions of the basic ones, this is not surprising and 
would occur in any quantum theory. 



4.2.4. Dynamics 
Also the constraints, in particular the Hamiltonian constraint 

written in Ashtekar variables with F^j^ = 2(3[„AJ^j + e'^*^AaA^, are non-polynomial functions of the basic variables and 
thus subject to quantization ambiguities. But given the difficulties of other attempts to formulate quantum gravity, even 
a single well-defined quantization of the Hamiltonian constraint would be a success. With the techniques described 
before, this can be accomplished 167^,6^. 

Also here, we need an inverse determinant of the triad which follows from the relation ( fTTI ). We have not encountered 
the curvature components before in operators, but they can be quantized by using the relation 

used also in lattice gauge theories. Here, a is a square loop of coordinate size A and with tangent vectors s", s" in 
a vertex as shown in Fig.[T3] One thus replaces the curvature components by holonomies around small loops which 
can then directly be represented on the Hilbert space. For this, one has to choose a prescription for the loops since the 
classical expression is simply evaluated in a point. The prescription gives rise to further quantization ambiguities which 
are, however, not as large as they would be for matter field theories on a background metric due to diffeomorphism 
invariance: only knotting and linking of the loop with edges in the graph of the state matters but not the precise position 
of an embedding in space. 

The remaining terms in the constraint involve extrinsic curvature and are the most complicated to deal with. 
We would have to subtract the spin connection from the basic Ashtekar connection to obtain extrinsic curvature 
components. Ashtekar connection components could again be quantized using holonomies, but the spin connection 
(fT2] l in general is a complicated function of the triad components. Moreover, it is not covariant, not even tensorial, and 
thus impossible to quantize directly. Instead of following this line it is conceptually easier to use a further identity II66II 




which allows one to express extrinsic curvature through commutators of the already quantized first term in the 
constraint with the volume operator and holonomies. The result is a highly complicated operator, but it is well-defined. 
Moreover, it displays crucial properties and deviations from the classical behavior which can be studied in models 
where the operator simplifies. 



4.2.5. Quantum effects 

The typical properties shown by the construction of quantum Hamiltonians are results of using holonomies, the 
basic tenet of loop quantum gravity: 

1. Quantized inverse powers of triad components give rise to modified small-scale behavior of coefficients. For 
singularities, this may be related to the issue of boundedness of coefficients in diff. equations as discussed before. 



2. Replacing local curvature and connection components by holonomies along extended loops implies non-locality 
as well as higher order spatial derivatives. This will be seen later to imply difference operators in equations for 
wave functions. 

In suitable semiclassical states, the quantum Hamiltonian must have an expectation value identical to the classical 
expression to zeroth order in /i. In any interacting (non-linear) theory, however, there will be quantum corrections which 
one can formulate in effective classical equations as will be discussed in the final section. In deep quantum regimes, 
full quantum equations have to be used with potentially very different properties compared to classical behavior or 
the Wheeler-DeWitt quantization as we will study them in the next section. This provides a general framework for 
quantum dynamics in which quantum hyperbolicity is testable. 

5. LOOP QUANTUM COSMOLOGY 

Loop quantum gravity provides a non-perturbative and background independent formulation of quantum gravity. Its 
main ingredients are a well-defined representation of basic fields, spatial discreteness and candidates for quantum 
dynamics. A description to study singularities from the perspective of quantum gravity is thus in principle provided 
but, in such a general setting, difficult to apply. There are not just severe complications from technical as well as 
conceptual issues of quantum gravity, but already the classical understanding of singularities in general is not precise 
enough even to decide where one would have to look for resolved singularities in quantum gravity. 

As in classical relativity, it is helpful in such a situation to study in detail explicitly treatable, usually symmetric 
situations where general aspects of quantum dynamics can be seen in action. This will at least give examples for 
singularity resolution and can suggest general mechanisms. By looking at different classes of models one then has a 
good chance of deciding whether or not such mechanisms are general or make use of special properties only realized 
in such models. Several unsuccessful examples are known, and also one so far successful scheme to be discussed here. 

When using symmetric models in quantum gravity one should be aware of differences between classical symmetric 
solutions, which are exact albeit special solutions of full general relativity, and symmetric quantum models. In contrast 
to classical solutions symmetric quantum solutions cannot be exact since uncertainty relations are violated. Both the 
configuration variables (the densitized triad E'^) and their momenta (the Ashtekar connection A[J must be symmetric 
to ensure a space-time solution respecting the symmetry everywhere. Non-symmetric modes of all canonical variables 
are thus zero which is possible classically but not in the presence of quantum uncertainty relations. A more general 
stability analysis is then required which has been done in a few cases. Alternatively, one can relax symmetries once 
models are well-understood and try to approach the full situation as closely as possible. 

An additional aspect of loop quantum gravity makes symmetric models worthwhile and shows that they can capture 
essential ingredients of a full quantization of gravity. Models can in particular illustrate consequences of quantum 
representations which have wide implications not only for the basic variables represented but for any composite 
operator constructed from them. The full representation induces distinguished representations of basic variables 
(analogously to holonomies and fluxes) in symmetric sectors which, at a kinematical level, are thus derived from 
the full theory |68]. Since the full representation is unique, the induced representations of models are distinguished 
among all representations one could try in a mini- or midisuperspace quantization. 

The induction proceeds by first identifying symmetric states in the full setting 1 69] . This is possible in a distributional 
sense, although the underlying discrete structure prevents the existence of any non-trivial normalizable state invariant 
under a continuous symmetry. The induced representation is then derived through basic full operators fixing these 
states, determining the reduced Hilbert space structure as well as the action of reduced basic operators (see the 
corresponding sections in |70,|2ll|72l). Composite operators then are to be constructed from those basic ones within 
the model following the same steps as in a full setting. A derivation from full composite operators would be more 
complicated and has so far not been attempted. But for testing properties of the basic loop representation the sketched 
procedure is sufficient. After all, dynamics has not been defined uniquely in the full theory and even an unambiguous 
derivation of the Hamiltonian constraint of a model from a full candidate would thus not result in a unique quantum 
dynamics. The situation in models is thus the same as that in the full theory where several candidates for dynamics 
formulated on a distinguished representation exist. The only difference is that models are often treatable explicitly and 
thus allow one to check physical consequences of different proposals for dynamics. In this way, the theory becomes 
physically testable. 

Models allow one to understand the full theory because its characteristic properties are preserved during the 
induction procedure. Most importantly, discreteness of spatial geometry is realized in an analogous way as we will 



see soon. Thus, even in isotropic models the representation is inequivalent to the Wheeler-DeWitt one which would 
have a continuous spectrum of the scale factor as a multiplication operator on the positive real line. This is the key 
reason why loop quantum cosmology, as the theory of symmetric sectors of loop quantum gravity is called, provides 
new insights even for the extensively studied field of quantum cosmology. 

In addition to the fact that quantum dynamics is often treatable explicitly in models, there is the added advantage 
of a much clearer classical singularity structure. In isotropic models, for instance, for the usual matter ingredients it 
suffices to formulate the condition a = to select singular states, which is easily done using the volume operator. 
Direct tests of quantum hyperbolicity then become possible. 



5.1. Isotropic models 

The basic quantities of an isotropic model formulated in Ashtekar variables comprise one conjugate pair {c,p) with 

\p\^y , c=i(^ + 7fl). (19) 

The only difference to metric variables is the fact that p can take both signs since it describes a densitized triad. The 
sign of p then is the intrinsic orientation of space. The parameter k in c is k — for a spatially flat model and k = I 
for positive spatial curvature. For negative spatial curvature the connection has a different form not covered here (see 

ill). 

The induced quantum representation must be a quantization of this finite dimensional model. However, it differs 
from the usual Schrodinger representation one would use in a Wheeler-DeWitt quantization. A complete set of 
orthonormal states is given by 1,75.1 

(c|;U) = e'>"/2 ^ (20) 

which akeady demonstrates the inequivalence to usual quantum mechanics where plane waves would not be normal- 
izable. Moreover, the Hilbert space is non-separable, i.e. has an uncountable basis. Basic operators, obtained by full 
flux and holonomy operators fixing the symmetric states ( l20b interpreted as distributions, act on these states by 

/3|M> - ^r^MlM) , ?^w5|m> = Im+m'>- (21) 

These operators indeed demonstrate the same properties we saw in the full theory: p has a discrete spectrum since 
all its eigenstates are normalizable, and e''"/^ is represented but not c itself which is the hallmark of any loop 
representation. (It is not possible to derive a c-operator by taking a derivative by ju' because the matrix elements of 
gi^i'c/i jjj j}jg basis states are not differentiable.) There is only a difference in how the discreteness of the spectrum of 
geometrical operators such as p is realized since the set of all eigenvalues of p is the whole real line. Nevertheless, 
the operator has a discrete spectrum since all its eigenstates are normalizable. In usual quantum mechanics, with a 
separable Hilbert space, this would imply that the set of eigenvalues is a discrete subset of the real line. On a non- 
separable Hilbert space as we have it here with an uncountable basis labeled by e M, however, an operator with a 
discrete spectrum can have any eigenvalue set. We will later see that the mathematical definition using normalizability 
of eigenstates is the one that is also relevant for the singularity issue. It is thus crucial that this property, rather than 
any statement about the set of eigenvalues, is preserved. This is a consequence of strong restrictions in the full theory 
from background independence and the subsequent transfer of the representation to symmetric models. 



5.7.7. Difference equation 

With the induced representation, we can follow most of the steps done in full to construct Hamiltonian constraint 
operators. The classical Hamiltonian constraint with contributions from a matter Hamiltonian 77matter is 

//= [7-Hc-k/2)^ + k^/4] y/\i\+Hramer{p,(p,P^)=0 (22) 

which can easily be seen to reduce to the Friedmann equation (with energy density p = Ipl^^/^T/matter) once c and p 
are replaced in terms of a and d using (fT9] l. This expression needs to be quantized by using "holonomies" e'^^'/^, with 
some /Xo G M to be chosen, for the connection components and Poisson brackets for the triad components. (Although 



the latter part is not necessary since there is no inverse of p due to cancellations from isotropy, we keep this step 
in order to have the quantization as close to the full one as possible.) The action of the resulting operatoil can be 
determined explicitly, choosing /Xo = 1 : 

{H-Hu..nM^j^^^-^{V^,+x-y^-l){e-''^^ (23) 

with the quantized matter Hamiltonian ^matter- In semiclassical regimes extrinsic curvature is small, c — k/2 <C 1, 
which by construction leads to the correct classical limit even though exponentials of c have been used. 

The operator equation i/| = to be solved for physical states can be expressed as a set of equations for expansion 
coefficients VV (^) I = L^i W (^) 1/^) which represent the state in the triad rather than connection representation. 
Applying the operator to such a general state and comparing coefficients of results in the difference equation 

- V^+3yVM+4(<P) - (2 + fc2)(y^+i -y^_i)rM('P) + - V^-5)e-'VM-4(<P) (24) 

= -57rG7'4^niatter(M)'/M('P) 

for \lffi{(p) written in terms of volume eigenvalues Vfi = {Y£p\^\/6)^^^ entering through commutators of holonomies 
with the volume operator 

This defines the dynamical equation for wave functions on minisuperspace, which can be appUed now especially 
in a neighborhood of the classical singularity at /i = in the interior. There are two sides to the classical singularity 
thanks to the triad orientation. The key question then is whether quantum propagation stops at = as the classical 
evolution would at p — 0. We have the matter Hamiltonian in a coefficient of the difference equation, which must be 
well-defined at = 0. This is our first test implied by quantum hyperbolicity in isotropic models as depicted in Fig.fTol 
The matter Hamiltonian Hfp = ^\p\^^^^p^ + \p\^^^V{(p), as before, contains p^^ which cannot be quantized directly 
since p has a discrete spectrum containing zero and thus lacks a densely defined inverse. 



5.7.2. Isotropic curvature bounds 

This is the place where the definition of a discrete spectrum through the normalizability of eigenstates is key since 
this, rather than properties of the set of eigenvalues other than zero being contained in the set, implies the non-existence 
of an inverse operator. It thus seems that the situation is even worse than in a Wheeler-DeWitt quantization since 
coefficients in the matter Hamiltonian appear impossible to define at all. But classically, a^^ can be rewritten in a form 
suitable for quantization, mimicking the full identity ( fTTI l of ll66ll . as 
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3/(2-2/) 



87r7G/;0-+l)(2; + l),t 



J^tijitMhjWpl'}) 



using only positive powers of p and "holonomies" hj — e"' . All this can directly be quantized, with the Poisson bracket 
becoming a commutator. Rewriting in this way introduces ambiguities because it can be done in many classically 
equivalent ways. Making different choices does, however, influence the quantization. Here, we have indicated two 
possibilities, taking different powers of V in the Poisson bracket or taking the trace in different representations. 
Eigenvalues of the resulting operators can be computed explicitly, 

/ ; \ 3/(2-2/) 

- — - (;■/)/ Q -' \ 

showing the dependence on the parameters j and I. 



There are different versions of this operator in the literature, which partially reflects the freedom existing in the full theory as well as details of 
the reduction to isotropic models which have not fully been worked out yet. For instance, one may reorder the operator, the one written here being 
non-symmetric, or include several effects which may be expected in an inhomogeneous quantization 1721 . 1761 . We will mention these possibilities 
here only when they are relevant for the singularity issue. 




Despite of the ambiguities, there are crucial common properties to all these quantizations. The divergence of the 
classical a^^ is cut off by quantum effects which is clear from the fact that the operator is well-defined on all basis 
states, even on |0) which corresponds to the classical singularity iIttIi . The expression for eigenvalues shows that 
inverse scale factor operators in fact all annihilate this state, irrespective of /' and I, since for = one is summing an 
odd expression in k from —j to j. One can most easily see the behavior around /i = by looking at an approximation 
of the eigenvalues valid for larger j. Viewing the sum as a Riemann sum of an integral, one obtains iItsIi 

with ii{p) =6p/Yil and 

P/(^) = |^^'-'(y^(('7+l)'+'-k-l|'+')-y^'?(('?+l)'+'-sgn(<?-l)|<7-ir+^)) . (26) 

Important properties, i.e. the approach to classical behavior at large jx, the peak at small values and the approach to 
zero for ^ = are all robust as can be seen from Fig. fl4ll.78[l79ll . 



5.1.3. Isotropic quantum hyperbolicity 

The difference equation ( l24l l for VV(^) msqA to investigate quantum hyperbolicity in isotropic models. We 

have already seen that the coefficients remain well-defined even at /i = 0. But this is only one first test which is 
necessary for quantum hyperbolicity to have a chance of being realized. We must, most importantly, be able to extend 
any wave function uniquely across = 0. For this, the recurrence scheme determined by the difference equation must 
be well-defined, i.e. we must be able to compute y/^ from the preceding values which themselves have been computed 
from some initial values at large /i. By following this procedure step by step one can see that evolution does continue 
from either side of the classical singularity (at /i > 0, say) to a new branch (at /X < 0) preceding the big bang at 
/X = 0. This scheme with two sides to a classical singularity is provided automatically by ingredients we were forced 
to assume in loop quantum gravity |80]. It is, in contrast to bounded curvature which was a condition only in isotropic 
models, a general scheme which is even realized in inhomogeneous models. Since we know the geometrical meaning 
of sgn(/i) which changes during the transition through the classical singularity, we can interpret the process as a 
change in orientation: the universe turns its inside out. By restricting variables to metrics one completely misses this 
possibility of non-singular behavior, but also the precise form of quantum dynamics is necessary for a well-defined 
transition. Around /X = 0, discreteness manifest in the difference equation is essential and classical space-time as a 
smooth manifold dissolves. 

The penetration of the classical singularity is thus non-trivial despite the discreteness: = is contained in the 
lattice and not simply jumped over. Moreover, leading coefficients of the difference equation may vanish, which could 
imply a break-down of the recurrence. In a backward evolution, for instance, we solve recursively for terms of 

y/^ and VV+4- This is only possible if V^-s — V|j_5, the coefficient of the difference equation, is nonzero. But 



it does vanish for jJ- — 4- since depends only on The value xj/o of the wave function at the classical singularity 
thus appears to remain undetermined. 

Nonetheless, xj/^ is determined uniquely for all positive and negative fi: i/o just decouples completely. We can follow 
the recurrence to negative values of fi. When determining yf-4, t/o seems necessary. But it drops out of the equation 
because then the coefficient Vi — V- 1 vanishes as well as the matter Hamiltonian which, as a robust property despite 
of quantization ambiguities, annihilates the state |0). The singular state |0) is then called manritQ with respect to the 
given evolution: It plays a passive role in the recurrence scheme of quantum evolution. In general, a mantic states can 
be defined as one at which the recurrence scheme implied by the Hamiltonian constraint in the triad representation 
changes its form. 

Mantic states have implications not just for the discussion of recurrence schemes, but they also imply dynamical 
initial conditions ll8Tl[82ll : Rather than determining y/Q which completely dropped out of the equations, the equation 
for ji = A implies a Unear relation between \//4 and \//g and thus an additional linear relation between initial values 
chosen at large /i > 0. There are thus restrictions on initial values not independent of but implied by the dynamical 
law, unlike the situation in every other area of physics. This implication, however, is more sensitive to the precise form 
of the constraint. Such conditions would be weaker for a symmetric ordering of the operator. 



5.2. Anisotropic quantum hyperbolicity 

As isotropic models are always very special, we have to drop symmetry conditions and see what remains of the 
observed mechanism of singularity removal. We first drop one of the isotropy conditions and look at space-times 
which are homogenous but have only one rotational symmetry axis. Such models are interesting for cosmology, but 
also for black hole physics since the interior inside the horizon of the Schwarzschild space-time is of this form. 

The densitized triad for the Schwarzschild interior can be written as 

in spherical coordinates where factors of sintJ arise due to the density weight of Ef . There are now two independent 
triad components, p^. and py. The determinant of the triad is det(£'") — Pcp\ and its orientation sgnp^ is solely 
determined by the sign of pc- The sign of ph is irrelevant, and in fact there is a residual gauge transformation pt, ~pb 
left after partially fixing the SU(2)-gauge by requiring the x-component of the triad to point in the su(2)-direction T3 
as used above. 

The triad determines a spatial metric 

\pA 

whose comparison with the interior Schwarzschild metric (r < 2m in the Schwarzschild coordinate where r becomes 
time-like and is called t in what follows), ds^ ~ {2m/ 1 — l)djic^ + f^di2^, allows us to identify the classical singularity 
on minisuperspace {p[,,pc): pc = at the Schwarzschild singularity while p;, = is not a singularity but the horizon. 

This model is loop quantized by a representation ph\ll,v) = \yt\ll\n,v), pc\n,v) = Y£pV\lJ.,v) of basic triad 
operators acting on orthonormal states |^,v) with ;U,v £ M, /X > 0. As in isotropic models, one can write the 
Hamiltonian constraint equation for states as a difference equation (using Xj/jj = ijf-^) llssll 

2 v/|v + 2| + Cf, ^\) v+2 - 2 (C^-2v/|v + 2| + v^) 1/^-2, v+2 
+ (v/|v+l|- v/|v-l|) ((Ai + 2)v/M+4.v-2(l +2f)^xj/^y + (Ai -2)v/^_4,v) 
+2(C^+2 v/|v-2| + v/^+2,v-2 - 2(C^-2 V\v-2\ + V/'^-2,v-2 = (27) 

where C^^\fi + 1/2| - 1/2| = sgn^. 

By the same procedure as before we conclude that this recurrence is singularity free which would also be realized 
if a matter term were present. Now, however, a more non-trivial test results: there are two boundaries of a metric 
minisuperspace, the horizon at = and the classical singularity at v = 0. Only one direction (v) is to be extended 



' "Manto: I stand still, around me circles time." Goethe's Faust 



in densitized triad variables if the singularity is removed. But there should be no extension through the boundary 
corresponding to the horizon because the validity of the homogeneous model breaks down outside the horizon. There 
is thus a non-trivial consistency check of the scheme: evolution through the classical singularity, but not the horizon 
must be realized. As one can see from the difference equation, the extension is provided in just the right manner By 
the general scheme, wave functions are extended only from one orientation of the triad to another, which provides the 
new branch at the other side of the classical singularity. Since orientation is determined by sgn{pc), only the boundary 
at = is penetrated, which corresponds to the classical singularity. The horizon pi, ~ 0, on the other hand, remains 
a boundary even for quantum evolution. That this is indeed non-trivial can be seen by trying to reproduce the scheme 
in co-triad or triad variables without the density weight. Although there is a new branch of the opposite orientation, 
the position of singularities and horizons in minisuperspace is different and no natural singularity resolution follows 



5.2.1. Beyond the singularity 

In general, solutions to difference equations, especially those of higher order as encountered in (l24l l and ( |27] i. could 
be wildly oscillating, such as xj/^ = (— 1)'' as a solution of VV+2 ~ + VV-2 — 0. Wave functions could thus be 
sensitive to Planck scale physics even in large volume or small curvature regimes where classical physics should be 
a good approximation. Even if initial values in a semiclassical regime are chosen to avoid this, oscillations could 
develop after evolving through a classical singularity. This would make it difficult to interpret the new branch as a 
classical one even when it extends to large volume. Moreover, oscillating solutions of difference equations can even 
be growing exponentially in amplitude and thus dominate any non-oscillating part. That this does not happen is a 
restrictive stability condition |85] which happens to be satisfied automatically in isotropic models as presented above. 

For the Schwarzschild interior, mantic states again play an important role for this issue. The coefficient i//2 o of the 
wave function at /i = 2, v = drops out of the highest order term, which implies the symmetry of solutions under 
V 1-^ —V [86]. If initial values at v > are chosen so as to suppress oscillations on small scales (giving a so-called 
pre-classical wave function), oscillations will not arise at the other side. This is in accordance with staticity we know 
is realized in the outside region of Schwarzschild: there should be no difference between the past and the future 
of the classical singularity. It is again highly non-trivial how this is realized here through mantic states around the 
classical singularity in minisuperspace, although only the interior is used which itself is not static. But the classical 
interior dynamics is determined by the same constraint as the outside, just restricted to the interior variables. Although 
the quantum dynamics used here is not derived from an inhomogeneous constraint but constructed within the model 
mimicking the construction of the full constraint, we arrive at quantum solutions in accordance with the classical 
expectation. The dynamics is responsible for the symmetry of solutions which was not required at the kinematical 
level, in agreement with the static outside as it is determined by the classical constraint. 

So far, it is unclear what happens in an inhomogeneous treatment of the black hole or with matter terms which 
would imply back-reaction on the geometry and non-trivial classical dynamics. A possible scenario concerning the 
information loss issue, based on several assumptions but taking into account singularity removal, has been described 
in iIstIi . Although suitable constraint operators and difference equations are available ISSll . they are very complicated 
to analyze. Some results are reviewed in [89, ,90] . 

Oscillations of solutions are also related to the normalizability of wave functions in an inner product. States not only 
have to solve the difference equation but, as in any quantum theory, must be normalizable. This is not always obvious 
to do when implementing a quantum constraint which usually changes the inner product on its solution space. In our 
discussion so far, and for most of the rest of this section, we can safely ignore this issue since we were able to show 
that all solutions are uniquely extended across classical singularities. This must then also be true for the physically 
normalizable ones. Here we see that it is important to show that all solutions, not just "generic" ones, are extended. 

Nevertheless, for more precise pictures of the evolution through a classical singularity knowledge of properties of 
the physical inner product can be relevant. One method to derive the inner product is group averaging [91[ which 
can be understood as writing quantum solutions to the constraint H\\lf) =0 as 8{H)\'^f) = f drexp(/fH)| y/). This is 
difficult to compute for the constraints we have to deal with here (especially due to the absolute value around jx which 
prevents the use of Fourier transformations, see e.g. [92^ for possible alternatives). But since H depends on momenta, 
its exponential is related to a shift operator in jj.. This suggests that normalizable states are indeed nearly constant 
on small scales and do not show strong oscillations. This justifies the condition of pre-classicality which has been 
introduced on intuitive physical grounds [81J. Oscillations only arise on larger scales where the matter term becomes 
important, or on smaller scales when curvature itself is large. Since, as we saw, pre-classicality can always be achieved 
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FIGURE 15. Effective curvature potential in a Bianclii IX model (47ll4^l48ll . While it is finite around the isotropy point in the 
center, despite of classical divergences in this region, it is unbounded when large anisotropies are reached. 



locally by choosing suitable initial conditions, but may be difficult to achieve on both sides of a classical singularity, 
strong restrictions on the quantization can be expected from a consistent physical normalizability. 



In anisotropic models, intrinsic curvature terms or the matter Hamiltonian occur in coefficients of the difference 
equation and must be well-defined everywhere for a consistent recurrence scheme. As seen before, in isotropy this 
implied finite inverse volume if quantum hyperbolicity is realized, but it is a very special case. Indeed we have a 
more general behavior in anisotropic models. Non-isotropic quantum hyperbolicity does not require boundedness of 
curvature in this sense, and it is indeed not realized generally as can be seen from Fig.fTSl 

But as we saw, the recurrence relations still do not break down at classical singularities and quantum hyperbolicity is 
realized. Bounded curvature on all of minisuperspace is thus not required for quantum hyperbolicity. In loop quantum 
gravity, unboundedness occurs in such a way that it does not prevent quantum hyperbolicity 1,51.1 . 



In inhomogeneous models we have not just one constraint equation but infinitely many ones since the Hamiltonian 
constraint has to be satisfied for any lapse function. Moreover, these are coupled equations although, for a wave 
function xj/, they remain linear For a spherically symmetric model, we have states 



with labels kg gWi for the edges and < /x,, G M for vertices. Again using a triad representation, coefficients of the 
wave function are subject to coupled difference equations (one for each edge) 



5.2.2. Unbounded curvature 



5.3. Inhomogeneous models 



(28) 



Co(^)V/(...,^-,^+,...)+C«+(^)V/(...,^-,^+-2,...)+C«_(fe)V/(...,fe-,fe++2,...) 
+Cl+(^)V/(. ..,k-- 2,k+, . . .)+CL-(k)xif{. . .,k^+2,k+, . . .) = 



(29) 



on an extended superspace of densitized triads. All coefficients Cj{k) are operators on the vertex labels jj. which we 
suppressed in the notation. They have all been computed explicitly in [88]. Local orientation sgndetZi is determined by 
sgnke such that we have to investigate the behavior of the coupled difference equations around vanishing edge labels. 



Again, evolution is non-singular |j42l which here depends crucially on the form (especially possible zeros) of 
coefficients Cr± (k) in a way which is much more non-trivial than in isotropic models. Unlike in homogeneous models, 
a symmetric ordering is required to extend solutions. Still the solution space is restricted by dynamical initial conditions 
as a consequence of mantic states. This shows that extending models to include more degrees of freedom does lead 
to tighter conditions on the allowed quantizations. So far one scheme is working in all situations considered, a highly 
non-trivial result given the complexity of equations such as ( |29] l compared to the much simpler case of ( |24] |. 

5.4. General properties 

The preceding examples exhaust all types of models where triad representations exist and which are so far treat- 
able explicitly. They allow possible general considerations: assuming singularities of BKL-type, only homogeneous 
behavior and diagonal metric or triad components are essential close to most interesting classical singularities. This 
would imply the existence of a triad representation at least for good approximations also around general singularities 
where the above arguments discussed in models would go through. That neither inhomogeneities nor local degrees 
of freedom by themselves spoil quantum hyperbolicity in loop quantum gravity follows from the demonstration that 
spherically symmetric and polarized cylindrically symmetric models respect the mechanism i42l 18811 . 

If this is not realized, non-commuting triad operators have to be taken into account. Metric components would thus 
be unsharp and the singularity appears washed out. An example realized in loop quantum gravity, as a consequence of 

non-Abelian SU(2)-holonomies, is discussed in ifsTll . While d{a) — a^^ used in (|25T l has the same eigenstates as d or 
the volume operator d^, this is not the case if full SU(2) holonomies rather than U(l) elements e''"^/^ are used. Using 
the basic equation (fTTI i. an inverse power is replaced by 

e"{Ai„V}^2e-hr{TMh;\v}) 

with a holonomy lig for some suitable edge e whose tangent vector is e" and whose parameter length is e. The right hand 
side is a good approximation for small e or small AJ,, which allows one to expand the exponential of the holonomy. 
The left hand side is not connection dependent, despite its appearance, since the derivative in the Poisson bracket 
removes A[,. For holonomies he taking values in an Abelian group, the right hand side is also connection independent 
since the two holonomies cancel each other even after taking a derivative by connection components. But due to 
non-commutativity, the cancellation is not complete for non-Abelian holonomies. Thus, the right hand side, which is 
used directly for a quantization, does become connection dependent. As a consequence, the resulting operator will not 
commute with the volume operator even though the classical expression does not depend on momenta of the densitized 
triad. 

This may lead to unsharp definitions of degeneracy points since it depends on which operator is used to determine 
eigenstates. There must be additional dynamical arguments to select the appropriate operator and its eigenstates as 
corresponding to classical singularities. Such situations would be much more complicated to analyze both from the 
classical perspective (structure of singularities) and for quantum dynamics. Fortunately, every indication so far, relying 
on dynamical properties of quantum gravity, points to properties as they were anticipated from models |51]. 

As we saw, symmetric models allow explicit investigations of many properties expected for singularities. It is 
important to keep a wide view of all types of characteristic models since any given model by itself may be too special as 
one can see it in isotropy. Still, there may always be properties not seen so far which may become relevant. According to 
the current status, quantum hyperbolicity is realized in much more general terms than any other mechanism to remove 
classical singularities. Moreover, it is suggested naturally by the structure of background independent quantum gravity 
(e.g. properties of singularities in densitized triads or the form of difference operators). Also the non-perturbative 
treatment matters, as shown in |70] where an anisotropic model was loop quantized not as described before but as 
a perturbation around a loop quantized isotropic model. Then, quantum hyperbolicity was not realized, in contrast 
to the full, non-perturbative quantization of the same model. Generalizations to less symmetric models have thus 
presented many non-trivial tests of the whole framework which so far were all passed successfully. A further test is 
the independence of the mechanism from details of the matter Hamiltonian, implying that we are dealing with a pure 
quantum geometry effect. Even curvature couplings, which can arise for non-minimally coupled scalar fields, do not 
change the mechanism |93] although at first the classical structure seems to be quite different from that in the absence 
of curvature couplings. The deep quantum picture crucially relies on spatial discreteness and dynamical equations for 
a wave function. This prevents detailed intuitive formulations, which in semiclassical regimes are sometimes available 



based on effective theory. Under certain conditions this can be applied to the transition through classical singularities, 
resulting in bounce pictures to be discussed in the next, final section. 

6. SEMICLASSICAL PICTURES 

As the general scheme of quantum hyperbolicity turns out to be realized in many cases in loop quantum gravity without 
counterexamples so far, it becomes possible to ask more detailed questions about the transition through classical 
singularities and in particular as to what happened before the big bang. The general mechanism utilizes in an essential 
way the discreteness of spatial geometry realized by the loop quantization: difference equations for wave functions on 
the space of triads need to be considered. 

Such equations are often difficult to analyze or to solve explicitly (see ll94l l95l l96l l97l losll for some techniques), 
and especially difficult to interpret. Moreover, at such a fundamental level, the issue of how space-time emerges from 
the underlying quantum state becomes exceedingly difficult especially on inhomogeneous models with many degrees 
of freedom. This is an issue faced by any quantum theory of gravity, independently of the precise methods used. For a 
recent discussion in string theory, using the AdS/CFT correspondence, see ll99ll . Intuitive pictures then usually require 
the use of special models allowing exact solutions, which imposes more than just a class of symmetries such as also 
specific matter ingredients (a free, massless scalar, say). While such models can often be analyzed in great detail, their 
properties may or may not be typical for general cases even within the same class of symmetries. This is analogous to 
the classical situation where the simplest and best known cases of singularities (isotropic cosmology, Schwarzschild 
black hole) may not display the general behavior of classical singularities as they are shown to arise from singularity 
theorems are from more general studies. Nevertheless, since these types of singularities are quite relevant for our 
understanding of the universe, detailed pictures for how they can be resolved are valuable. Moreover, they can provide 
a basis for more general scenarios, for instance in a perturbation theory. The most common intuitive picture for such 
resolutions is a bounce 



6.1. Effective equations 

Much intuition in cosmology derives from isotropic or at least homogeneous models where, from a semiclassical 
perspective, a bounce is the only way to resolve a classical singularity. There are only finitely many (at most three) 
metric components in those models which must stay away from zero to avoid a singularity. Thus, within a finite 
amount of time all components evolve through a minimum if they were initially contracting and thus lead to a bounce 
in volume. Note that this argument makes use of the finite dimensionality of the classical phase space and does thus not 
apply to inhomogeneous models. Moreover, the argument presupposes that some kind of semiclassical description at 
least of the quantum gravity state (if not matter fields) is valid throughout, which can then be identified with a smooth 
spatial geometry subject to an altered, bouncing dynamics. 

While these assumptions are difficult to justify in general, or to be avoided in bounce arguments, bounce models can 
be useful to study semiclassical effects of quantum gravity. Moreover, if sufficiently well developed, they can be used 
to understand how cosmological inhomogeneities could evolve through a bounce with potential effects on structure 
formation. Rather surprisingly, one particular bounce model provides the zeroth order basis ("free theory") for an 
effective theory of loop quantum gravity and a corresponding perturbation scheme. The technical details are fully 
analogous to the widely used low energy effective actions in particle physics and thus very promising for developing 
a detailed understanding of the semiclassical properties of loop quantum gravity. 

6.7.7. Low energy effective action 

It is often most useful to describe quantum effects in certain regimes by correction terms to classical equations, 
without changing the type of differential equations (except for possible higher derivative terms). Thus, a quantum 
mechanical system would still be formulated in terms of coupled ordinary differential equations in time rather than a 



Some exa mples for bounces from loop effects can be found in llO(lllOllll02l[T03lll04 . l74h . Oscillatory scenarios have been developed and 
described in fTo a n M flOT i n M fTol. 



partial differential equation. Or, a quantum field theory would be described by partial differential equations in space- 
time coordinates rather than functional differential equations. This does not only lead to technical simplifications in 
solution procedures but also alleviates most interpretational issues of quantum theories. 

In quantum mechanics, it is well-known that all harmonic oscillator states follow classical trajectories exactly. There 
is thus no need for modifications to classical equations of motion, and it turns out that only a zero point energy is added 
to the Hamiltonian Hq = jj^p^ + jin^ax^q^ + jhco. This is similar for free quantum field theories although the zero- 
point energy may diverge. For an-harmonic contributions U to the potential or for interactions, there are corrections 
to the classical Hamiltonian which can be computed in a perturbation theory around the harmonic oscillator or a free 
theory. By the usual particle physics techniques based on Legendre transforms of generating functionals of n-point 
functions this leads to the low-energy effective action 
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showing first order corrections in h to the kinetic as well as potential term of the classical action. 

These low energy effective actions as they result from perturbations around a vacuum state need to be generalized 
for quantum gravity. In that case, the Hamiltonian is not only unbounded from below but also does not allow a general 
notion of energy to supply meaning to "low energy effective action." (For quantizations of gravity on a background 
manifold one can and often does employ low energy effective actions. While this is valuable for scattering effects, it 
is not suitable and can be highly misleading in the context of cosmology.) Such generalizations are available and can 
be derived, for instance, from a geometrical formulation of quantum mechanics: One can view the Hilbert space of a 
quantum system as an infinite-dimensional phase space whose Poisson brackets are given through the imaginary part 
of the inner product. This satisfies all mathematical conditions imposed in the definition of Poisson brackets. 

To derive the Poisson structure in a simple way, we choose expansion coefficients in states = LyCyl'^/) with 
respect to an arbitrary basis { | xj/j) } as coordinates on the Hilbert space. Then, the imaginary part of the inner product 
of states can be interpreted as Poisson structure 

{Rec^-,Imci.} = ^^Sjk 

equipping the Hilbert space with a phase space structure. Furthermore, the Schrodinger equation for | y/) turns out to 
be equivalent to Hamiltonian equations of motion generated by the quantum Hamiltonian //g({cy}) = (H). Here, Hq 
is a function on the phase space with coordinates cj by using the state with expansion coefficients cj to compute the 
expectation value of H. That this gives the correct dynamics follows easily by choosing the eigenbasis { | } of H. 



Then,//Q({c;}) =l,jc,c;( 
equations of motion are 



\j/i\H\\j/j) — with energy eigenvalues Ej, H\\jfj) — Ej\\jfj), and the Hamiltonian 



d E 
-Rscj = {Rscj^Hq} = ylmc; 



d E 
-Imcj = {lmcj,HQ} = -yRecy . 



(31) 



This is equivalent to cj = —ih^^Ejcj in agreement with the Schrodinger evolution. Quantum mechanics is thus 
formally much closer to classical mechanics, making a relation through effective equations possible in suitable 
approximations. While the full quantum dynamics already appears in classical form for the infinitely many variables 
Cj, any effective approximation must include a truncation to finitely many variables for a mechanical system. 



6.7.2. Quantum variables 

A more useful set of coordinates for this purpose is defined as follows: we use "classical" variable^ 

q = {q) and p = {p) 



(32) 



This term often gives rise to confusion as these variables are close to the classical ones (e.g. in the sense of describing the peak position of a 
single wave packet) only in semiclassical regimes while they are defined here for an ai'bitrary quantum system. The justification for the term is that 
irresjj ective of the regime the quantum phase space is a fiber bundle over the classical phase space with bundle projection 1 1;/) i//), i//)) 



and quantum variables 

G"-":=((^-(^)r«(/3-(/3)r)weyi (33) 

where n > 2, a = 0, . . . ,« and "Weyl" denotes the totally symmetric ordering of the operators before the expectation 
value is taken. Poisson relations of these variables are related to commutators as we will see in more detail later: 
q,p} = { \q, p\)/ih = 1 and {q,G" "} = ~ {p,G"'"}. For {G"'",G^''"} a closed formula exists but is rather lengthy 

imitni. 

Quantum variables are dynamical just as the classical variables are: for a semiclassical state they change in time, e.g. 
if a wave packet spreads and deforms. The resulting evolution back-reacts on the classical variables which determine 
the peak position of the wave packet. Their motion then in general differs from the classical one, which is to be captured 
in appropriate quantum correction terms to the classical equations of motion. The exact behavior is determined by the 
Schrodinger equation, or equivalently by the quantum Hamiltonian Hq which couples classical and quantum variables. 
For instance, for a cubic potential we have (q^) =q^ + 6qG^-^ + 60^'^ with a coupling term qG^'^ in addition to a zero- 
point contribution by G"'^. 

More generally, for an an-harmonic oscillator with classical Hamiltonian H = 4^^p^ + jtnco^q^ + U [q] we have a 
quantum Hamiltonian with infinitely many coupling terms, 

Hq = {H{q,p)) = {H{q+iq-q),p + {p-p))) (34) 

2m 2 2 'ToW' \m(0 I 

written in dimensionless variables &'•" = fi^"^^{mco)"^^^"G"'" . As indicated, these terms follow from formally 
expanding the Hamiltonian in ^ — ^ and p — p- Using the Poisson brackets of all classical and quantum variables, 
Hq generates Hamiltonian equations of motion / = {/,//g}: 

q ^ - (36) 



m 



p = -,nC0\-U'{q)-j:U^) U^"+'\q)&- (37) 



G"'" = -a(0&'-^'" + in~a)(0&'+^"-a--^&'-^'" (38) 
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These are infinitely many coupled equations for infinitely many variables. At this stage, the system is still fully 
equivalent to the Schrodinger dynamics, just written in terms of "n-point functions" q, p and G" " instead of the 
wave function they determine. 

As an example we can look at the harmonic oscillator whose Hamiltonian equations of motion are 

p= {p,Hq} =-map-q 



1 

q= {q,HQ\ 



m 

G«'«= {G"-",//„} = -(«-fl)G"+''"-m«2flG"-''". 
m 

In this case, all terms coupling the G" " for different n vanish, and we have an infinite set of differential equations 
only finitely many of which are coupled to each other Moreover, the equations are linear and can easily be solved. 
For instance, constant solutions for the quantum variables exist, satisfying uncertainty relations such as G^'^G^'^ > 

^ + (G ' '2 ) 2 . These solutions correspond to the well-known coherent states which do not spread and follow the classical 
trajectories exactly. But even if quantum variables are not constant, they do not appear in equations of motion for the 
classical variables and thus do not back-react on them. This is why the effective action of the harmonic oscillator is 
identical to its classical action. 



With an-harmonic contributions, coupling terms are switched on and all equations get coupled to each other Consis- 
tent truncations to finitely many equations for finitely many variables are then required for an effective approximation. 
This is possible in, e.g., an adiabatic approximation ill 1 ill : we solve approximately for the leading G"" assuming that 
they change in time much more slowly than q and p, and insert the solutions into the equations of motion for q and p. 
Doing this to first order in h and to second order in the adiabatic approximation, and writing the first order equations 
for q and p as a second order equation for q, we obtain 



nu"'{q) 
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+mOp-q + U'{q) + ^-^ ryr = 

as it also follows from the low energy effective action (jSOjl. 



6.1.3. Solvable systems 

The derivation clearly shows the role of the harmonic oscillator: its classical and quantum variables decouple to 
sets of finitely many linear equations and there is no back-reaction by quantum variables on expectation values. This 
happens whenever the Hamiltonian is quadratic in canonical variables, or more generally when the system is linear, 
i.e. when the Hamiltonian and a set of basic variables form a Lie algebra using Poisson brackets. Effective equations 
can then be obtained by perturbing around the exact solutions of such a solvable system. 

Realizing this opens one possibility for generalizations of low energy effective equations, which are then based on 
alternative solvable systems suitable for a given context. Remarkably, examples for such solvable models are realized 
in cosmology and allow a systematic effective theory in this context. 



6.2. Large scale effective theory for cosmological bounces 

Identifying such a solvable model suitable for cosmology and deriving its perturbation equations allows one to 
derive intuitive semiclassical pictures which describe the transition through classical singularities in detail. 



6.2.7. Solvable model for cosmology 

The harmonic oscillator with its periodic motion is not suitable as a solvable system for cosmology. But we can 
look at a free isotropic scalar model whose Friedmann equation — ^^p^^^^p^ follows from (l22l) with variables 

c = a (extrinsic curvature) and p^^^ = a^ (volume) where p is the densitized triad component, now assumed positive. 
Solving for p(p yields p,p oc cp —: H, to be interpreted as the Hamiltonian which generates the flow in the variable (p 
playing the role of internal time. This Hamiltonian is quadratic|3 although not of the harmonic oscillator form. But as 
in this case, classical and quantum variables decouple and the quantum Hamiltonian Hq = cp + is obtained, as in 
( |35] |, by adding only the zero point contribution G'^'' ~ j{{cp) + (pc) ) — cp. (Here, we use a slightly different notation 
compared to the general quantum variables ( |33] | for better clarity.) 

Its equations of motion follow by using Poisson brackets {c,G'''} = = {/9,G''^} as well as {G" , G'^} — IG'^'^, 
{G" , G'^'^l = AG'^P, {G^'^jG''''} — 2G'''' and further ones depending on which variables one is interested in solving 



We were not careful about the signs involved when solving the quadratic constraint equations for p^, although one can conclude from the 
constraint only that \ « \cp\ . Using the absolute value for H, on the other hand, would not leave it strictly quadratic. As we are mainly interested 
in states for which the expectation value of H is large compared to the spread AH, we do not need to worry about significant contributions from 
solutions with different signs of the Hamiltonian. Note that H and AH are preserved in time. Thus, if the condition H S> AH is satisfied once, e.g. 
for an initial semiclassical state, it will be satisfied at all values of (p. 



for. We will now show how the relevant Poisson relations are derived, using the example of {G^'^jG''^}. The basic 
identity is the relation 

{{A),{B)} = {[A,B])/ih (39) 

which clearly shows how the quantum Poisson brackets are related to commutators. The familiar relation {q,p} = 
{[q,p]) /iTi = 1 then follows immediately. But it cannot be directly applied to quantum variables such as G" — (c^) — 
and G''' = ^{{cp) + {pc))—cp since here also products of expectation values occur This can easily be dealt with using 
the Leibniz rule 

{f,8i82}^8i{f,82} + {f,8i}82 (40) 

to reduce all Poisson brackets to those of expectation values of operators where ( [39] ) applies. With the brackets 

{(c^), (c/3)} = 2(c2), {{c^),cp} = {(c2), (c)(p)} = c{(c2), (/3)} = 2c^ and {c^,cp} = 2c^ we then derive 

{G"-,G''^} = - c2, i((c/3) + (/3c)) - cp} = 2((c2) - c') = 2G- . 

Similarly, all other Poisson brackets are derived, for which also closed formulas exist 111 llll I2I1 . 

Equations of motion generated by Hq = cp + C^'' are thus c — c,p — —p for the classical variables, and — 2G'^'^, 
G^P — and G'''' — —2G'''' for the quantum variables. This is easily solved by c{t) — cje', p{t) — C2e^' , G"'{t) = cge^', 
G'^P{t) —C4 and G''''{t) —c^e^^' with suitable integration constants which are only restricted by the uncertainty relation 
C3C5 > tP'/4 + c^. Although constant solutions of the quantum variables do not exist, the semiclassical properties 
are quite similar to those of the harmonic oscillator In particular, semiclassicality is preserved: the relative spreads 
^p/p = V GPP I p and Ac/c = \fG^ jc are constant throughout the whole evolution. This will allow us to derive 
effective equations and to develop a perturbation theory around the known solutions of the solvable model where the 
free isotropic scalar plays the same role as the harmonic oscillator in quantum mechanics or free theories in quantum 
field theory. It also shows that conclusions drawn from the free isotropic scalar have to be taken with great care, just as 
the behavior of the harmonic oscillator is not at all typical for general quantum systems even in semiclassical regimes. 



6.2.2. Loop formulation and bounces 

In a loop formulation we do not use the Schrodinger quantization of basic variables c and p. Instead, in a 
loop quantization (ISTT i the operator p has a discrete spectrum and no operator for c exists. What is represented 
are only exponentials exp(/c) such that, e.g., sine occurs instead of c. The resulting Hamiltonian operator H = 
— j!(exp(/c) — exp(— /c))/3 is a shift operator and implies a difference equation for the state in a triad representation. 
This operator is not identical to (|23] ) derived earlier but is closely related. What we have not included here are quantum 
effects in the inverse power p^^^^ of the matter Hamiltonian as they occur in ( |25] |. Such terms do not allow solvable 
models but can be included in perturbation theory. 

The Hamiltonian is non-quadratic and not solvable in an obvious way. But introducing J — pe"^ allows us to reorder 
the Hamiltonian to become a linear expression H = —ji{f—P). The price to pay is that the algebra of basic operators 
p and J is non-canonical, which usually implies that the system is not solvable in the above sense even for a linear 
Hamiltonian. But for the syst em u nder consideration it turns out that the set of Hamiltonian operator and basic variables 
{p,J) forms a linear system ll 13ll . given by the (trivially) centrally extended sl(2,R) algebra 

\p,f\=hj , [p,P] = -riP , [f,P]^ -2hp~h^ . (41) 

Taking expectation values, the linear quantum Hamiltonian Hq = —ji{J — J) does not even receive a zero point 
contribution and generates equations of motion 

p = {p,HQ}^-^iJ + j) , J^{J,HQ}^-{ip + h)^j (42) 

with solution 

p{t)^^{cie-'+C2e')~^h , J{t) = ^{cie-'-C2e') + iH. (43) 

They are simply linear combinations of the solutions in the Schrodinger quantization. Depending on whether ciC2 is 
positive or negative we have bouncing solutions of cosh form or solutions of sinh form which arrive at the classical 
singularity p — Q after a finite amount of internal time t. For our purpose of studying singularity removal we have to 
analyze these types of solutions in more detail. 
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FIGURE 16. Bouncing effective solutions for expectation value and spread. 



But we have not yet implemented all necessary conditions and some of those solutions are non-physical. Classically 
we have JJ = for J — pe\p{ic), which as a reality condition for c must have an analog in the quantization. For 
states, this corresponds to normalizability in a physical inner product ensuring that exp(/c) is quantized to a unitary 
operator Although computing a physical inner product is usually a difficult issue, for the solutions p{t) and J{t) of 
expectation values we can implement the reality condition directly by just noticing that JJ — p^ + 0{h) must remain 
satisfied at the effective level. Only corrections of order h to the classical condition may result since (J) (P) ^ {Jf')- 
The only way to implement this at all times is by requiring 

JJ ^ {cie-'-C2e'f+4H^ ^ 
{p + Ti/2)^ {cxe-'+C2e'Y ^ ^ 

which implies ciC2 = up to quantum coiTections which are not important for large H. This leaves only the bouncing 
solution 

p{t)=Hcosh{t-5)-h , y(f) = -//(sinh(f-5)+/) (45) 

with a single constant of integration 8. Note that the minimum of p is given hy H — h which for a large Hamiltonian 
H, i.e. large matter content, is far away from the classical singularity at p — 0. The bounce trajectory agree s well wit h 
numerical solutions of (physically normalized) wave packets in a closely related model studied recently in lll04lll4l . 
Similarly, we can compute equations of motion for the spread parameters 

GPP = -2G"^ , G-'^ = -2G"^ , ^-^G-'-' -^GPP-^ip'--Jj+np + h^/2) (46) 
(using Jp — p^ and the commutation relations). They satisfy the uncertainty relation 

qPPqJJ _^QpJ^2^_^j^2^ 

For H ^ h, a solution is given by (A/?)^ = G^f « fiH cosh{2{t — §2)). The semiclassical behavior throughout the 
bounce is clearly seen from Fig. [16] where during the contraction and expansion phases the relative spreads are almost 
constant as in the Schrodinger quantization, although they may change from the contracting to the expanding phase 
depending on the integration constant 5 — 52. The loop quantization thus connects the two branches in a well-defined 
way but leaves open the relative degrees of semiclassicality on both sides. 



6.2.3. Perturbations 

As described in general terms before, perturbations around solvable models can be formulated by standard means. 
This allows one to introduce terms which imply couplings between classical and quantum variables, such as matter 
potentials, inhomogeneous degrees of freedom, or basic issues of a loop quantization. The latter effects include 



different factor orderings of the constraint or modifications due to inverse powers in the matter Hamiltonian. This 
provides a systematic way to derive the evolution of fields through cosmological bounces, which is not available in 
other schemes where assumptions about the regularity at a bounce must be imposed rather than being derivable. The 
perturbation scheme mentioned here can certainly break down, which implies that the semiclassical nature of realistic 
bounces is testable in a self-consistent manner. 

In addition to conceptual lessons, mechanisms for singularity resolution provide a possib le rel ation to observations 
in bounce scenarios of structure formation, making use of Hamiltonian perturbation theory ill 15ll and an implementa- 
tion in perturbative loop quantum gravity 111 16(l . 

7. THE CURRENT STATUS OF SINGULARITY RESOLUTION IN QUANTUM 

GRAVITY 

We have focused the discussions here on the singularity problem from a general perspective. Accordingly, detailed 
scenarios we described were only those which have been formulated in a sufficiently general context, applying to 
more than a small class of situations and not showing any obvious limiting assumptions. 

With recent developments, loop quantum cosmology has provided the first systematic effective theory for the 
non-singular evolution of perturbations through a classical singularity using perturbation theory around bouncing 
solutions. Expectation values and the spreading of semiclassical states can be computed explicitly. This shows under 
which circumstances semiclassical states are obtained also during and after the bounce, although the degree of 
semiclassicality depends on initial conditions. A general analysis is unfinished. In general, bounce pictures are difficult 
to extend to inhomogeneities unless the latter remain perturbative. Often, one appeals to the BKL conjecture to suggest 
that homogeneous scenarios might be sufficient to discuss singularity resolution. But by construction bounces in 
general do not allow one to go sufficiently deep in the regime, asymptotically close to a classical singularity, which is 
used in the BKL scenario to argue that time derivatives dominate spatial derivatives. 

Under unrestricted perturbations, bounce pictures can easily break down and full quantum properties are required. 
Even in this case, singularity removal in loop quantum gravity is established in many situations by quantum hyperbol- 
icity: wave functions on spaces of triads can be extended uniquely across classical singularities. This crucially rests 
on quantum geometry (rather than matter) and especially its discrete spatial structure. Background independence is 
important for the detailed realization, providing many consistency tests for quantum gravity. Quantum hyperbolicity 
is at present the most general mechanism for non-singular behavior, indeed showing clearly the role of important as- 
pects of a background independent quantization of gravity. Nevertheless, much remains to be understood for a general 
statement on singularity resolution. The main problem at the current stage of developments is not only the complexity 
of quantum gravity, but also an understanding of the classical structure of singularities needed to select, in the absence 
of general solutions, the quantum regime to be looked at. 

ACKNOWLEDGMENTS 

The author is grateful to Mario Novello for an invitation to give a lecture course at the Xllth Brazilian School on 
Cosmology and Gravitation on which this article is based. He also profitted from talks and discussions presented at the 
workshop "The Quantum Nature of Spacetime Singularities" held at the KITP, Santa Barbara, in January 2007, and 
thanks Claus Kiefer for comments on possible singularity removal in a Wheeler-DeWitt quantization. This research 
was supported in part by the National Science Foundation under Grant Nos. PHY99-07949 and PHY05-54771. 

REFERENCES 

1. M. Bojowald,|http : / / relativity . phys . Isu . edu/ilqgs /bo jowald022 007 . pdf| 

2. R. M. Wald, General Relativity, The University of Chicago Press, 1984. 

3. H. Friedrich, and A. D. Rendall, Lect. Notes Phys. 540, 127-224 (2000),'gr-q c/0002 074 l 

4. S. W. Hawking, and R. Penrose, Pwc. Roy Soc. Land. A 314, 529-548 (1970). 

5. S. W. Hawking, and G. F. R. Ellis, The Large Scale Structure of Space-Time, Cambridge University Press, 1973. 

6. J. M. M. Senovilla, Gen. Rel. Grav. 29, 701-848 (1997). 

7. R. Penrose, Phys. Rev. Lett. 14, 57-59 (1965). 



8. B. Krishnan, and E. Schnetter, Phys. Rev. D 73, 021502(R) (2006), .gr-qc/0511017| 

9. I. Ben-Dov, igr- qc / 0611057 

10. J. M. M. Senovilla, Phys. Rev. Lett. 64, 2219-2221 (1990). 

11. F. J. Chinea, L. Fernandez- Jambrina, and J. M. M. Senovilla, Ph ys. Rev. D 4 5, 481-486 (1992), [gr^qc/ 04 03075] 

12. C. J. S. Clarke, Class. Quant Grav. 15, 975-984 (1998), gr-qc/ 97 02 033j 

13. J. A. Vickers, and J. P. Wilson, Class. Quant. Grav. 17, 1333-1360 (2000), gr-qc/ 9907 105| 

14. J. A. Vickers, and J. P. Wilson, gr-qc/0101018 

15. A. Horde, A. H. Guth, and A. Vilenkin, Phys. Rev. Lett. 90, 151301 (2003), [gF^qc / 0110012] 

16. M. R Djbrowski, Annalen Phys. (Leipzig) 15, 352-363 (2006). 

17. G. T. Horowitz, and R. C. Myers, Gen. Rel. Grav 27, 915-919 (1995), jgr-qc/ 95030 62] 

18. J. Polchinski, String Theory, vol. I and II, Cambridge University Press, 1998. 

19. B. Craps, hep-th/0605199 

20. A. Strominger, Nucl. Phys. B 451, 96-108 (1995), [hiF-th/9504 0"90] 

21. G. T. Horowitz, and A. R. Steif, Phys. Rev. Lett 64, 260-26 3 (1990). 

22. H. Liu, G. Moore, and N. Seiberg, JHEP 0206, 045 (2002), l hep-th/0204168 

23. H. Liu, G. Moore, and N. Seiberg, JHEP 0210, 031 (2002), hep -th/0206182 

24. J. McGreevy, and E. Silverstein, JHEP 0508, 090 (2005), hep-thT oSO 6130| 

25. R. Brandenberger, and C. Vafa, Nucl. Phys. B 316, 391 (1989). 

26. M. Bojowald, and A. Skirzewski, "Quantum Gravity and Higher Curvature Actions," in Current Mathematical Topics in 
Gravitation and Cosmology ( 42nd Karpacz Winter School of Theoretical Physics), International Journal of Geometric 



Methods in Modem Physics, 2006, hep-th/0606232 




27. G. T. Horowitz, and J. Polchinski, Phys. Rev D 66, 1035 


12 (2002), hep-th/02 622 8 


28. G. Festuccia, and H. Liu, JHEP 0604, 044 (2006), "hep- 


th/050620^ 



29. T. Hertog, and G. T. Horowitz, JHEP 0504, 005 (2005), he p-th/0503071| 

30. D. L. Wiltshire, 'An introduction to quantum cosmology," in Cosmology: The Physics of the Universe, edited by B. Robson, 
N. Visvanathan, and W. S. Woolcock, World Scientific, Singapore, 1996, pp. 473-531, gr- qc/0101003 

31. B . S . DeWitt, Phys. .Rev. 160, 1 1 1 3- 1 1 48 ( 1 967) . 

32. V. A. Belinskii, I. M. Khalatnikov, and E. M. Lifschitz, Adv Phys. 13, 639-667 (1982). 

33. G. F R. Elhs, and M. A. H. MacCallum, Commun. Math. Phys. 12, 108-141 (1969). 

34. D. Garfinkle, Phys. Rev Lett. 93, 161101 (2004), gr-qc/0312117 

35. A. D. Rendall, "The nature of spacetime singularities," in 100 Years of Relativity - Space-Time Structure: Einstein and 
Beyond, edited by A. Ashtekar, World Scientific, Singapore, 2005, igr-qc/0503112| 

36. J. B. Hartle, and S. W. Hawking, Phys. Rev D 28, 2960-2975 (1983). 

37. A. Vilenkin, Phys. Rev D 30, 509-5 1 1 (1984). 

38. H. D. Conradi, and H. D. Zeh, Phys. Lett. A 154, 321-326 (1991). 

39. P Hajicek, Nucl. Phys. B 603, 555-577 (2001), hep-th/ 0"007 005| 

40. P Hajicek, Led. Notes Phys. 631, 255-299 (2003), gr-qc/02 04049| 

41. M. Bojowald, Class. Quantum Grav 20, 2595-2615 (2003), gr- qc/0303073| 

42. M. Bojowald, Phys. Rev Lett. 95, 061301 (2005), gr-qc/0506 128| 

43. M. Bojowald, Class. Quantum Grav 18, L109-L1 16 (2001), gr -qc/0105113| 

44. M. Bojowald, and H. A. Morales-Tecotl, "Cosmological applications of loop quantum gravity," in Proceedings of the Fifth 
Mexican School (DGFM): The Early Universe and Observational Cosmology, Lect. Notes Phys. 646, Springer- Verlag, 2004, 
pp. 421-462, gr-qc/0306008 

45. R. Gambini, and J. Pullin, Phys. Rev Lett. 90, 021301 (2003), [gr^qc/02 6055| 

46. R. Gambini, and J. Pullin, Class. Quant Grav 20, 3341-3364 (2003), gr-qc/ 02 12 03 31 

47. M. Bojowald, G. Date, and K. Vandersloot, Class. Quantum Grav 21, 1253-1278 (2004), gr-qc/0311004| 

48. M. Bojowald, and G. Date, Phys. Rev Lett 92, 07 1302 (2004), gr-qc / 0311003 

49. M. Bojowald, G. Date, and G. M. Hossain, Class. Quantum Grav 21, 3541-3569 (200 4), gr-qc/0404039] 

50. A. Ashtekar, A. Corichi, and J. Zapata, Class. Quantum Grav 15, 2955-2972 (1998), gr-qc/9806041| 

51. M. Bojowald, Class. Quantum Grav 23, 987-1008 (2006), gr-qc/0508118 

52. R. Arnowitt, S. Deser, and C. W. Misner, "The Dynamics of General Relativity," in Gravitation: An Introduction to Current 
Research, edited by L. Witten, Wiley, New York, 1962. 

53. K. V. Kuchaf, Phys. Rev D 4, 955-986 (1971). 

54. B. Berger, Phys. Rev D 11, 2770-2780 (1975). 

55. A. Ashtekar, Phys. Rev D 36, 1587-1602 (1987). 

56. J. F Barbero G., Phys. Rev D 51, 5507-5510 (1995), gr- qc/ 941 0014| 

57. G. Immirzi, Class. Quantum Grav 14, L177-L181 (1997). 

58. R. Giles, Phys. Rev D 24, 2160-2168 (1981). 

59. C.Rovem,andL.Smolin,/'V-^ev. £152,5743-5759 (1995). 

60. J. Lewandowski, A. Okolow, H. Sahlmann, and T. Thieman, Commun. Math. Phys. 267, 703-733 (2006), [gr^qc/ 504147] 

61. C. Fleischhack, math-ph/0407005 

62. C. Fleischhack, Phys. Rev Lett. 97, 061302 (2006). 

63. C. Rovelh, and L. Smolin, Nucl. Phys. B 442, 593-619 (1995), erratum: Nucl. Phys. B 456, 753 (1995), [gF^qc/ 94 11005] 



64. A. Ashtekar, and J. Lewandowski, Class. Quantum Grav. 14, A55-A82 (1997),[ gr-qc/ 9602 04 6| 

65. A. Ashtekar, and J. Lewandowski, Adv. Theor. Math. Phys. 1, 388^29 (1997), g r-qc/9711031| 

66. T. Thiemann, Class. Quantum Grav. 15, 839-873 (1998), gr-qc / 960 608 9. 

67. C. Rovelli, and L. Smolin, Phys. Rev. Lett. 72, 446^49 (1994), gr-qc/93o'8002| 

68. M. Bojowald, Living Rev. Relativity 8, 11 (2005), 

|http: //relativity .livingreviews ■org/Articles/lrr-2005-l_l/|[gr-qc/0 6_01085| 

69. M. Bojowald, and H. A. Kastrup, Class. Quantum Grav. 17,"30CI9-3043l2000), 'hep-th799'6"7 Qg? 

70. M. Bojowald, H. H. Hernandez, and H. A. Morales-Tecotl, Class. Quantum Grav 23, 3491-3516 (2006), [gr^qc / 0511058] 

71. M. Bojowald, Class. Quantum Grav 21, 3733-3753 (2004), gr-qc/0407017 

72. M. Bojowald, Gen. Rel. Grav 38, 1771-1795 (2006), gr-q c/0609034 

73. M. Bojowald, Class. Quantum Grav 19, 2717-2741 (2002), |gr-qc / 02 02 77] 

74. K. Vandersloot, gr-qc/ 612 07 

75. A. Ashtekar, M. Bojowald, and J. Lewandowski, Adv Theor Math. Phys. 7, 233-268 (2003), gr -qc/0304074| 

76. A. Ashtekar, T. Pawlowski, and P. Singh, Phys. Rev D 74, 084003 ( 2006), igr-qc/0 607 03 9| 

77. M. Bojowald, Phys. Rev D 64, 084018 (2001), gr-qc 7oT"050 67| 

78. M. Bojowald, Class. Quantum Grav 19, 5113-5130 (2002), gr-qc/0206053' 

79. M. Bojowald, "Loop Quantum Cosmology: Recent Progress," in Proceedings of the International Conference on Gravitation 
and Cosmology (ICGC 2004), Cochin, India, Pramana, 2004, vol. 63, pp. 765-776, |gr-qc/04 02053| 

80. M. Bojowald, Phys. Rev Lett. 86, 5227-5230 (2001), gr-qc/0102069 

81. M. Bojowald, Phys. Rev Lett. 87, 121301 (2001), gr-qc/0104072 

82. M. Bojowald, Gen. Rel. Grav 35, 1877-1883 (2003), gr-q c/ 0305 0"69] 

83. A. Ashtekar, and M. Bojowald, Class. Quantum Grav 23, 3 9T-4 11 (2006), |g r - qc/0509075| 

84. L. Modesto, Class. Quantum Grav 23, 5587-5601 (2006), gr-qc/05090 78| 

85. M. Bojowald, and G. Date, Class. Quantum Grav 21, 121-143 (2004), gr- qc/0307083| 

86. D. Cartin, and G. Khanna, Phys. Rev D 73, 104009 (2006), gr-qc / 0602 025^, 

87. A. Ashtekar, and M. Bojowald, Class. Quantum Grav 22, 3349-3362 (2005), gr-qc/ 050402 9| 

88. M. Bojowald, and R. Swiderski, Class. Quantum Grav 23, 2129-2154 (2006), gr-qc/0511108 

89. M. Bojowald, "Quantum Riemannian Geometry and Black Holes," Trends in Quantum Gravity Research, Nova Science, 
2006, gr-qc/0602100 

90. M. Bojowald, Int. J. Mod. Phys. D 15, 1545-1559 (2006), proceedings of "Einstein's Legacy in the New Millenium," Puri, 
India, December 2005, gr -qc/0607130 

91. D. Marolf, Refined algebraic quantization: Systems with a single constraint, |gr-qc/ 9508015| 

92. B. Bahr, and T. Thiemann, gr-qc/0607075| 

93. M. Bojowald, and M. Kagan, Class. Quantum Grav 23, 4983^990 (2006), gr-qc/0604105| 

94. D. Cartin, G. Khanna, and M. Bojowald, Class. Quantum Grav 21, 4495-4509 (2004), gr -qc/0405126| 

95. M. Bojowald, and A. Rej, Class. Quantum Grav 22, 3399-3420 (2005), gr-qc/050 4100| 

96. D. Cartin, and G. Khanna, Phys. Rev Lett. 94, 111302 (2005), gr-qc/ 050 1 1 6 

97. G. Date, P/jy^. i?ev. 72, 067301 (2005), gr-qc/050503^ 

98. J. Rosen, J.-H. Jung, and G. Khanna, Class. Quantum Grav 23, 7075-7084 (2006), [gr^qc/0 607 044] 

99. G. Festuccia and H. Liu, h ep-th/0 6 110 98 

100. P Singh, and A. Toporensky, Phys. Rev D 69, 104008 (2004), 'j F^qc / 3 1 2 1 l"o] 

101. M. Bojowald, R. Maartens, and P Singh, Phys. Rev D 70, 08351 7 (2004),|hep-th/0407115l 

102. G. Date, and G. M. Hossain, Phys. Rev Let t 94, 011302 (2005 ), |gr-qc/0407074| 

103. G. Date, Phys. Rev D 71, 127502 (2005), gr-qc/050500_2| 

104. A. Ashtekar, T. Pawlowski, and P Singh, Phys. Rev Lett. 96, 141301 (2006), gr-qc/0602086| 

105. J. E. Lidsey, D. J. Mulryne, N. J. Nunes, and R. Tavakol, Phys. Rev D 70, 063521 (2004), gr-qc704 "0 6042[ 

106. D. J. Mulryne, N. J. Nunes, R. Tavakol, and J. Lidsey, Int J. Mod. Phys. A 20, 2347-2357 (2005), gr-qc/0411125J 

107. D. J. Mulryne, R. Tavakol, J. E. Lidsey, and G. F. R. Elhs, Phys. Rev D 71, 123512 (2005), astro-ph/0502589. 

108. M. Bojowald, Nature 436, 920-921 (2005). 

109. N. J. Nunes, Phys. Rev D 72, 103510 (2005),'astro-ph/0507 683l 

1 10. A. Ashtekar, and T. A. Schilling, "Geometrical Formulation of Quantum Mechanics," in On Einstein's Path: Essays in Honor 
of Engelbert Schiicking, edited by A. Harvey, Springer, New York, 1999, p p. 23-65,igr-qc / 97 60 6 9| 

111. M. Bojowald, and A. Skirzewski, i?ev. Math. Phys. 18, 713-745 (2006), math-ph/0511043| 

112. A. Skirzewski, PhD Thesis, HU Berlin, 2006. 

113. M. Bojowald, gr-qc / 0608100 

1 14. A. Ashtekar, T. Pawlowski, and P Singh, Phys. Rev D 73, 124038 (2006), gF^7 6 6 OTT 

115. M. Bojowald, H. Hernandez, M. Kagan, P Singh, and A. Skirzewski, Phys. Rev D 74, 123512 (2006),|gr-qc/0609057| 

1 16. M. Bojowald, H. Hernandez, M. Kagan, and A. Skirzewski, Phys. Rev. D to appear, g r-qc/ 0611112] 



